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13.  ABSTRACT  (Usxlmum  200  words) 

In  the  context  of  design  studies  for  the  U.S.  Large  Blast/Thermal  Simulator  (LB/TS),  the  TJ.  S. 
Army  Ballistic  Research  Laboratory  (BRL)  has  sponsored  optical  studies  of  the  flow  start-up  in 
convergent-divergent  nozzles  which  have  the  flow-initiating  diaphragm  located  in  the  nozzle  throat. 
The  experiments  were  performed  in  the  200-mm  shock  tube  at  the  Ernst-Mach-Institutc  in  Freiburg 
(Breisgau),  West-Gcrmany.  The  scope  of  the  studies  included  divergent  nozzles  with  ha’f-cone  angles 
of  6°,  16°,  45°  and  90°  tested  at  seven  diaphragm  pressure  ratios  ranging  from  four  to  188  and  result¬ 
ing  in  shock  strengths  ranging  from  1.4  to  4.4.  The  results  were  summarized  in  graphs  of  significant 
parameters,  e.g.,  shock  formation  time,  flow  start-up  period,  flow  expansion  angle  and  shock  strength 
versus  the  driver  pressure  ratio  and  further  compared  with  numerical  results  obtained  with  the  BRL- 
Q1D  hydrocode.  The  numerical-experimental  comparison  shows  qualitative  agreement.  The  flow 
phenomena  are  generally  reproduced  by  the  computations  as  long  as  they  are  not  strongly  dependent 
on  viscous  effects.  The  study  shows  that  a  pressure  loss  of  10  percent  is  connected  to  the  presence  of  a 
large  area  discontinuity  at  the  exit  plane  of  the  nozzle  throat  when  no  divergent  nozzle  is  attached. 
The  results  suggest  that  a  45°  divergent  nozzle  may  present  an  acceptable  compromise  for  minimizing 
these  pressure  losses  by  reducing  the  associated  area  discontinuities. 
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1.  INTRODUCTION 


It  is  known  that  the  gas  behind  a  shock  travelling  down  a  shock  tunnel  can  be 
expanded  to  high  flow  Mach  numbers  by  passing  it  through  convergent-divergent 
(DeLaval)  nozzles.  Such  nozzles  have  been  extensively  investigated,  both  in  theory 
and  experiment,  and  a  summary  of  these  efforts  has  been  presented  by  Amann 
(1968).  However,  in  all  of  these  studies  the  diaphragm  which  separates  the  high- 
pressure  region  from  the  low-pressure  region  was  located  upstream  of  the  nozzle 
under  investigation  so  that  a  well-formed  shock  entered  the  nozzle  from  the 
upstieam  side.  No  efforts  were  undertaken  to  study  the  shock  formation  in  nozzles 
immediately  after  flow  initialization. 

Large  Blast/Thermal  Simulators  (LB/TS)  are  specialized  shock  tunnels  for 
simulating  blast  and  thermal  effects  associated  with  decaying  blast  waves  such  as  are 
generated  by  nuclear  explosions.  LB/TS’s  are  postulated  to  employ  convergent- 
divergent  nozzles  at  the  driver  exit  to  retard  the  outflow  of  the  high-pressure  driver 
gas  thus  generating  long  flow  durations  (Cadet  and  Monzac  1981;  Mark,  Opalka,  et. 
al.  1983;  Pearson,  Opalka  and  Hisley  1985;  Opalka  and  Pearson  1988,  1989).  The 
nozzle  design  for  an  LB/TS  differs  from  the  nozzle  design  used  in  the  previously 
mentioned  studies  in  that  it  has  the  flow-initiating  diaphragm  mounted  in  the  throat 
of  the  nozzle  itself.  When  the  diaphragm  is  ruptured,  a  shock  forms  travelling 
downstream  into  the  divergent  section,  while  a  rarefaction  fan  travels  upstream  into 
the  convergent  section.  No  previous  experimental  research  is  known  for  this  prob¬ 
lem. 

The  U.S.  Army  Ballistic  Research  Laboratory  (BRL),  in  cooperation  with  the 
Ernst-Mach  Institut  (EMI),  Freiburg  im  Breisgau,  Germany,  therefore,  initiated  an 
optical  study  of  the  flow  start-up  process  in  four  convergent-divergent  nozzles  which 
had  the  diaphragm  separating  the  high-pressure  and  low-pressure  sections  of  the 
shock  tube  mounted  in  the  throat  of  the  nozzles.  Under  the  auspices  of  a  bi-lateral 
data  exchange  agreement  between  the  governments  of  the  U.S. A.  and  the  Federal 
Republic  of  Germany  (FRG),  the  experiments  were  performed  in  the  200mm  shock 
tube  at  EMI.  The  diaphragm  was  located  at  the  upstream  edge  of  the  viewing  win¬ 
dow  of  high-quality  optical  glass  through  which  shadow  and  schlieren  photographs 
could  be  taken.  The  optical  measurements  were  supplemented  by  measurements  of 
the  static  pressure  at  three  locations  downstream  from  the  nozzles. 

The  objectives  of  this  investigation  were  to  (a)  obtain  optical  records  of  the 
shock  formation  in  the  convergent-divergent  nozzle  immediately  after  rupturing  the 
diaphragm,  (b)  study  the  shock  formation  and  flow  start-up  processes  in  the  diver-  j 

gent  nozzles,  (c)  determine  the  influence  of  the  cone  angle  and  length  of  a  divergent  j 

nozzle  on  the  pressure  signature  in  the  expansion  region  behind  the  nozzle  and  (d) 
facilitate  comparisons  with  one-  and  two-dimensional  hydrocode  computations.  The 
test  setup  and  the  results  of  this  experimental  study  are  presented  in  this  report  and 
the  significant  findings  are  discussed.  ] 
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2.  BACKGROUND 


The  U.S.  Army  has  proposed  the  construction  of  a  test  facility  for  simulating 
ideal  blast  waves  and  large  enough  to  test  full-sized  military  equipment  in  order  to 
meet  its  growing  need  for  blast  and  thermal  survivability  testing.  The  BRL  was 
assigned  the  lead  role  in  the  research  and  development  phase  of  this  project.  Subse¬ 
quently,  the  BRL  developed  a  concept  of  such  a  facility,  referred  to  as  LB/TS,  over 
the  past  nine  years.  In  1988,  the  Defense  Nuclear  Agency  (DNA)  assumed  respon¬ 
sibility  for  the  construction  of  such  a  facility  to  be  located  on  the  White  Sands  Mis¬ 
sile  Range  (WSMR),  New  Mexico  (Figure  1).  The  proposed  U.  S.  LB/TS  is  suitable 
to  simulate  both  thermal  and  blast  effects  of  nuclear  explosions  over  a  wide  variety 
of  shock  overpressures  and  weapon  yields  without  generating  nuclear  radiation 
effects  (Figure  2).  The  LB/TS  will  be  used  to  conduct  blast  and  thermal  survivabil¬ 
ity  testing  on  full-scale  military  equipment,  like  trucks  and  helicopters,  and  to  per¬ 
form  research  into  nuclear  blast  phenomenology. 

2.1  LB/TS  Concept.  An  LB/TS  is  basically  a  shock  tunnel  whose  cross-sectional 
area  varies  along  its  length.  Constant  area  (i.e.,  straight)  shock  tubes  do  not  pro¬ 
duce  the  flow  durations  which  are  observed  in  decaying  blast  waves.  To  obtain  long 
flow  durations  comparable  to  those  of  decaying  blast  waves,  the  outflow  of  the 
driver  gas  has  to  be  retarded.  This  is  effected  by  converging  the  flow  area  at  the 
exit  of  the  driver,  so  that  the  driver  lakes  the  shape  of  a  bottle.  Furthermore,  a  sin¬ 
gle  driver  is  technically  impractical  because  of  the  size  of  the  facility  and  the 
required  supply  pressures  in  the  driver.  A  number  of  smaller  drivers  has  to  be 
employed  to  accomodate  the  necessary  initial  test  conditions.  The  BRL  LB/TS  con¬ 
cept  is  described  in  detail  by  Pearson  (1987)  and  Opalka  (1987,  1989)  and  illustrated 
in  Figure  3. 

Blast  waves  will  be  simulated  by  releasing  compressed  gas  from  several  high- 
pressure  steel  driver  tubes  into  a  large  expansion  tunnel  constructed  of  prestressed 
concrete.  The  gas  exits  the  drivers  through  convergent  nozzles.  Fast-acting  valves 
were  to  be  mounted  in  the  throats  of  these  nozzles;  but  DNA  abandoned  the  use  of 
valves  for  flow  initiation  in  favor  of  the  proven  diaphragm  technique.  The  diver¬ 
gent  nozzles  were  eliminated  for  cost  efficiency.  The  expansion  tunnel  has  a  semi¬ 
circular  cross  section  with  a  10  m  radius.  This  size  is  deemed  necessary  to  avoid 
blocking  of  the  flow  about  the  largest  target.  The  test  section  for  'he  targets  is  cen¬ 
tered  at  about  seven  hydraulic  diameters  (14.94  m)  downstream  from  the  exit  of  the 
driver  nozzles  in  the  expansion  tunnel.  The  thermal  simulation  will  be  effected 
through  aluminum/oxygen  combustion  near  the  target.  An  active  reflection-wave 
eliminator  (RWE)  at  the  open  end  of  the  expansion  tunnel  prevents  reflected  shocks 
or  rarefaction  waves  from  travelling  back  into  the  test  section  and  disfiguring  the 
shape  of  the  simulated  blast  wave. 

2.2  LB/TS  Design  Studies.  The  BRL-Q1D  code  (Opalka  and  Mark  1986)  was  used 
to  perform  computational  studies  for  determining  the  design  parameters  of  the  pro¬ 
jected  facility.  Parametric  studies  (Opalka  1987,  1989)  were  initiated  to  answer 
questions  about  the  necessary  size  and  the  expected  performance  of  such  a  facility 
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for  simulating  the  required  range  of  blast  waves.  The  scope  of  the  computational 
studies  encompassed  pressure  ratios  from  4  to  200  across  the  diaphragm  for 
unheated  driver  gas  with  stagnation  temperature  T4  =  288. 15  K  and  ratio  of  specific 
heats  y=  1.4  in  association  with  the  operating  envelope  for  the  planned  U.S. 
LB/TS  (Figure  2).  Standard  atmospheric  conditions  were  assumed  as  ambient  in  the 
expansion  tunnel  for  all  computational  simulations.  Table  I  presents  the  test  condi¬ 
tions  which  were  selected  for  the  parametric  BRL-Q1D  computations.  The  shock 
strength,  the  positive-phase  duration,  and  the  static  and  dynamic  impulses  of  the 
simulated  blast  wave  in  the  test  section  of  the  LB/TS  were  defined  from  the  initial 
test  conditions,  i.e.,  the  driver  volume,  V4,  the  driver  pressure  ratio,  P4/P1,  and  the 
driver  temperature  ratio,  T4/Tj,  where  subscript  1  refers  to  driven  (i.e.,  ambient) 
and  subscript  4  to  driver  conditions.  For  the  purpose  of  this  study,  the  initial  throat 
shock  overpressure,  pr^,  was  determined  by  iteratively  solving  the  basic  shock  tube 
equation,  P4/P1  =  f(PT2/Pi),  where  subscript  2  identifies  conditions  behind  the 
incident  shock. 


TABLE  I.  LB/TS  INITIAL  TEST  CONDITIONS 

P!  =  101.325  kPa;  Tt  =  288.15  K 

Pi,  psi 

Pi ,  kPa 

2.0 

13.79 

5.0 

34.47 

TEST  SECTION 

10.0  15.0  20.0 

68.95  103.4  137.9 

25.0 

172.4 

30.0 

206.8 

35.0 

241.3 

P4/P1 

P4 ,  MPa 

4 

0.4053 

14 

1.419 

DRIVER 

38  55  80 

3.850  5.573  8.106 

110 

11.145 

145 

14.692 

188 

19.050 

PT2/P1 

P7i ,  kPa 

1.929 

94.11 

3.256 

228.6 

THROAT 

4.689  5.288  5.974 
373.8  434.5  504.4 

6.575 

565.0 

7.120 

620.1 

7.652 

674.0 

Initially,  the  U.S.  design  studies  were  based  on  the  Large  Blast  Simulator 
(LBS)  at  the  Centre  d’Etudes  de  Gramat  (CEG),  France  (Gratias  and  Monzac 
1981).  However,  it  soon  became  clear  that  a  much  larger  facility  was  needed  in 
order  to  accommodate  the  full  range  of  anticipated  targets  (Ethridge,  et.  al.  1984). 
Also,  a  broader  range  of  shock  overpressures  and  weapon  yields  was  required  to 
cover  the  test  conditions  specified  in  the  operating  envelope  (Figure  2).  The  results 
of  our  computational  studies  with  the  BRL-Q1D  code  showed  that  the  full-scale 
blast/thermal  simulation  facility  must  employ  driver  gas  heating,  an  active  RWE, 
and  should  use  computer-controlled  throat  valves  to  optimize  control  and  operation 
of  the  facility.  The  divergent  nozzles  were  eliminated  from  the  design;  however, 
some  questions  remain  about  the  usefulness  of  such  divergent  nozzles. 

Experimental  efforts  were  initiated  at  the  same  time,  but  test  results  became 
available  only  much  later  because  the  construction  of  the  model  facility  required  a 
long  time.  To  date,  small-scale  experiments  with  heated  driver  gas,  as  well  as  an 
active  RWE  model,  have  been  completed  and  the  results  have  been  used  to  validate 
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the  computational  predictions  (Coulter  1987;  Gion  1989;  Schraml  and  Pearson 
1990).  A  throat-valve  model  has  been  built  and  is  awaiting  testing  for  the  validation 
of  the  design  concept  and  the  computational  predictions.  Also,  a  1/6-scale  test  bed 
using  a  double  diaphragm  and  an  active  RWE  is  in  the  construction  stage  and,  when 
built,  will  be  used  for  validating  the  design  concepts,  for  research,  and  for  vulnera¬ 
bility  testing  of  small  items  of  military  equipment. 

2.3  Nozzle  Flow  Phenomena.  The  drivers  of  the  proposed  U.  S.  LB/TS,  as  well  as 
of  the  CEG-LBS,  have  the  diaphragm  which  separates  the  high-pressure  region  from 
the  low-pressure  region  of  the  shock  tunnel  mounted  in  the  throat  of  a  convergent 
nozzle.  The  CEG-LBS  uses  divergent  nozzles  to  control  the  expansion  of  the  flow 
into  the  expansion  tunnel.  When  the  diaphragm  is  ruptured,  a  shock  travels  down¬ 
stream  into  the  divergent  section,  while  a  rarefaction  fan  travels  upstream  into  the 
convergent  section.  The  time  from  the  rupturing  of  the  diaphragm  to  the  establish¬ 
ment  of  a  quasi-steady  flow  pattern  in  the  nozzle  is  defined  as  the  starting  process  of 
the  nozzle.  No  previous  experimental  research  is  known  for  either  the  starting  pro¬ 
cess  or  the  quasi-steady  flow  phase  of  these  nozzles. 

The  divergent  nozzles  were  eliminated  from  the  U.  S.  design  because  they  seem 
to  present  certain  problems.  Figure  4  presents  an  experimental  and  a  computational 
pressure  history  of  a  42  kPa/10  kT  blast  wave.  It  is  argued  that  the  pressure  spike 
that  appears  at  the  front  of  the  blast  wave  in  the  experimental,  as  well  as  in  the  com¬ 
putational  pressure  history,  is  caused  by  the  divergent  nozzles.  This  argument  is 
supported  by  computational  results  (Opalka  1987)  which  indicate  that  the  greater  the 
cone  half-angle  of  the  divergent  nozzle,  the  smaller  is  the  spike  at  the  shock  front. 
It  is  argued  also  that  the  divergent  nozzles  increase  the  thrust  on  the  reaction  pier 
because  the  driver-nozzle  combination  functions  like  a  rocket  motor. 

Eliminating  the  divergent  nozzles  apparently  seems  to  have  a  double  advantage 
on  the  design  and  the  performance  of  the  U.S.  LB/TS  reducing  the  thrust  on  the 
pier  as  the  drivers  empty,  as  well  as  eliminating  the  spike  at  the  front  of  the  simu¬ 
lated  blast  wave.  On  the  other  hand,  results  of  small-scale  pilot  experiments  (Gion 
1989)  seem  to  indicate  that  there  are  pressure  losses  connected  to  the  introduction  of 
a  large  area  discontinuity  at  the  exit  of  the  nozzle  throat  which  are  lessened  when 
divergent  nozzles  are  used.  Therefore,  the  use  of  divergent  nozzles  cannot  be  easily 
dismissed.  However,  if  divergent  nozzles  are  included  in  the  U.S.  LB/TS  design, 
the  nozzle  length  becomes  a  critical  material  parameter  which  itself  depends  of  the 
cone  half-angle  chosen  for  the  nozzle  design. 

Amann  (1968)  has  shown  in  his  experiments  that  the  start-up  time  is  shortest 
when  the  cone  half-angle  of  the  divergent  nozzle  is  largest.  However,  his  research 
was  limited  to  cone  half-angles  of  5°-15°,  and  it  is  not  known  whether  or  not  the 
trend  continues  at  higher  angles.  The  BRL  one-dimensional,  and  two-dimensional 
flow  computations  (Opalka  1987;  Hisley  and  Molvik  1988)  a^ce  with  earlier 
theoretical  and  experimental  research  that  a  recompression  shock  forms  in  the  noz¬ 
zles  and  either  stays  there  for  moderate  diaphragm  pressure  ratios,  or  is  swept  down 
the  expansion  tunnel  for  higher  diaphragm  pressure  ratios.  For  want  of  knowledge. 
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and  for  comparison  with  the  numerical  predictions,  experimental  information  is 
needed  on  the  flow  start-up  process  with  and  without  divergent  nozzles  from  15°  to 
90°,  on  the  shock  formation  period,  the  formation  and  location  of  recompression 
shocks,  and  on  the  effect  of  the  cone  half-angle  on  the  quasi-steady  flow  once  the 
start-up  process  is  completed. 


3.  EXPERIMENTAL  APPARATUS 

The  experimental  apparatus  at  the  EMI  (Reichenbach  and  Opalka  1990)  con¬ 
sists  of:  A  shock  tube,  which  was  especially  modified  for  use  in  this  study;  a  24- 
spark  Cranz-Schardin  camera  and  shadowgraph  arrangement;  and  electronic  control 
and  recording  equipment. 

3.1  EMI  Shock  Tube.  The  EMI  shock  tube  has  a  circular  cross  section  with  an 
inner  diameter  of  200mm.  The  driver  section  is  1.800m  long,  and  the  driven  sec¬ 
tion  measures  6.980m  in  length  (Figure  5).  It  is  followed  by  a  1.570m  long  exten¬ 
sion  of  rectangular  crossection,  which  houses  the  test  chamber.  The  rectangular 
cross  section  of  the  test  chamber  is  110mm  high  by  40mm  wide  with  its  center 
located  7.210m  downstream  from  the  diaphragm.  The  side  walls  of  the  test 
chamber  are  formed  by  plane  parallel  windows  (200mm  long  by  110mm  high)  of 
high  quality  optical  glass  to  facilitate  optical  investigations  of  the  shock  wave  propa¬ 
gation.  An  inlet  channel  of  the  same  shape  and  cross  section  as  the  test  chamber 
and  formed  of  polished  aluminum  plates  protrudes  641mm  upstream  into  the 
axisymmetric  driven  section.  A  rectangular  section  is  cut  out  of  the  arriving  shock 
wave  by  the  sharp  inlet  edges  of  the  channel  walls  (cooky  cutter  principle)  to 
prevent  the  boundary  layers  which  have  formed  inside  the  driven  section  from 
entering  the  test  area. 

The  existing  shock  tube  had  to  be  modified  in  order  to  carry  out  the  special 
task  requested  with  this  project.  The  modifications  include  a  new  driver  pressure 
chamber  which  can  be  mounted  inside  the  inlet  channel,  and  the  various  model  parts 
for  the  divergent  nozzles  (Figure  6).  The  length  of  the  driver  pressure  chamber  was 
chosen  to  be  1  m  in  order  to  delay  the  arrival  of  reflections  from  the  end  wall.  The 
new  driver  pressure  chamber  includes  the  convergent  section  and  that  part  of  the 
nozzle  throat  which  lies  upstream  of  the  diaphragm  (see  Figure  6a).  The  inner  cross 
section  of  the  driver  is  64mm  high  by  30mm  wide;  it  converges  to  a  cross  section 
that  is  16mm  high  by  30mm  wide  at  the  diaphragm.  Note  that  the  diaphragm  open¬ 
ing  shown  in  Figure  6b  is  larger  than  the  hight  of  the  nozzle  throat.  It  is  reduced  to 
the  hight  of  the  throat  by  lips  on  the  half-nozzle  models  fitting  into  the  opening  (see 
Figure  7).  The  diaphragms  are  composed  of  two  Ultraphan  sheets,  each  0.1mm 
thick,  with  an  ignition  wire  of  0.08  mm  diameter  sandwiched  between  them  (Figure 
6b).  Cellulose  acetate  sheet  (Ultrapna'.)  has  the  desirable  qualities  of  becoming 
brittle  when  stretched  and  breaking  rapidly.  With  this  construction,  the  diaphragm 
can  be  burst  at  an  exactly  predetermined  moment  in  time  by  exploding  the  embed¬ 
ded  wire  with  a  high-voltage  pulse. 
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The  test  objectives  require  that  the  diaphragm  be  mounted  in  or  near  the  test 
section  window  of  the  shock  tube  in  such  a  way  that  the  view  of  the  downstream 
region  of  the  throat  and  divergent  nozzle  is  only  minimally  obstructed.  The  lengths 
of  the  various  divergent  nozzle  configurations  (see  Figures  6a  and  7)  necessitate 
relocating  the  models  in  the  test  section  window  so  that  optical  records  may  be  taken 
of  all  waves  of  interest  from  one  (200mm)  to  three  window  lengths  (600mm) 
downstream  from  the  diaphragm.  This  requirement  makes  it  necessary  to  locate  the 
diaphragm  either  at,  or  one  or  two  times  the  window  length  upstream  from,  the  left 
edge  of  the  test  section  window.  The  window  locations  are  labeled  B,  C,  and  D  in 
Figure  6a.  The  location  A  was  not  used  in  this  experimental  effort. 

3.2  Nozzle  Configurations.  Four  convergent-divergent  nozzle  models  were  fabri¬ 
cated  (Figure  7),  thus,  producing  the  four  different  layouts  of  the  divergent  sections 
shown  in  Figure  6a.  All  dimensions  were  scaled  proportionately  to  the  area  ratios 
of  the  convergent-aivergent  nozzles  used  in  the  CEG-LBS.  The  45°  convergent  sec¬ 
tion  is  identical  for  all  nozzle  models  and  part  of  the  driver  assembly.  The  four 
divergent  nozzle  configurations  are  characterized  as  follows: 

a)  A  90°  divergent  nozzle,  i.e.,  a  throat  without  a  divergent  section  dumping  the  gas 
directly  into  the  expansion  tunnel  (Insert  2  shown  in  Figures  6a  and  7a). 

b)  A  divergent  nozzle  with  a  45°  cone  half-angle  extended  to  the  full  height  of  the 
expansion  tunnel  (Inserts  2  and  3  shown  in  Figures  6a  and  7b). 

c)  A  16°  divergent  nozzle  with  the  area  and  length  ratios  of  the  Q  ID-computational 
CEG-LBS  model  (i.e.,  the  7  CEG-LBS  nozzles  combined  into  one)  extended  to  the 
full  height  of  the  expansion  tunnel  (Inserts  2  and  4  shown  in  Figures  6a  and  7c). 

d)  A  6°  divergent  nozzle  with  the  area  and  length  ratios  and  the  cone  half-angle  of 
the  nozzles  used  in  the  CEG-LBS  design  (Inserts  2  and  5  shown  in  Figures  6a  and 
7d). 

3.3  Instrumentation.  A  24-frame  Cranz-Schardin  spark  camera  and  shadow- 
schlieren  system  was  the  principal  tool  for  the  optical  study  (Figure  8).  The  sparks 
from  24  point  sources  are  focussed  onto  24  objective  lenses  of  the  camera  by  means 
of  a  concave  mirror.  At  the  moment  a  spark  is  generated,  it  is  projected  onto  the 
film  plane  through  the  objective  lens  onto  which  the  image  of  the  spark  is  copied. 
This  procedure  has  the  advantage  of  needing  no  moving  parts  to  align  the  path  of 
the  light  rays.  Therefore,  the  optical  resolution  is  determined  essentially  by  the 
aperture  of  the  objective  lenses.  Parallax  may  present  a  disadvantage  in  some  cases. 
However,  the  effect  of  parallax  is  minimal  in  the  present  setup  where  a  concave  mir¬ 
ror  with  a  long  focal  length  (f  =  3.5m)  is  used  in  conjunction  with  a  test  section  of 
shallow  depth  (40mm).  Parallax  was  found  in  these  tests  to  be  less  than  one  percent 
of  the  length  scale. 

The  density  gradients  in  the  flow  are  made  visible  by  the  schlieren  method. 
The  basic  idea  of  the  schlieren  method  is  that  part  of  the  deflected  light  is  inter¬ 
cepted  before  it  reaches  the  photographic  plane,  so  that  the  areas  of  the  flow  field 
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through  which  it  has  passed  appear  darker.  The  Toepler  technique  (Schardin  1934), 
using  sharp  knife  edges,  was  found  to  yield  the  best  results  (see  Figure  8).  The  sha¬ 
dow  method  visualizes  the  second  derivative  of  the  density  and  is  well  suited  for  the 
visualization  of  shock  waves.  The  knife  edge  is  removed  and,  to  generate  the  sha¬ 
dow  effect,  the  objective  lenses  are  focussed  on  a  reference  plane  typically  at  a  dis¬ 
tance  of  400  mm  beyond  the  center  plane  instead  of  on  the  center  plane  of  the  test 
chamber.  A  point  source  of  light  is  required  for  this  method  to  work  well.  In  the 
EMI  system,  the  light  source  has  a  diameter  of  less  than  1mm.  In  the  following 
experiments,  both  methods  were  combined  and,  in  this  way,  clear  images  were 
obtained  of  the  shock  waves  as  well  as  of  the  density  gradients. 

Other  equipment  used  included  a  high-voltage  power  unit,  the  trigger  gauges 
and  electronic  delay  units  for  the  timing  of  the  Cranz-Schardin  camera,  and  several 
12-bit  transient  recorders.  Kistler  quartz  pressure  gages  of  type  603B  were  used  to 
sense  the  static  pressure.  They  have  a  characteristic  frequency  response  of  approxi¬ 
mately  500  kHz.  Electrical  filters  with  a  maximum  frequency  response  of  180  kHz 
were  installed  between  the  pressure  gauges  and  the  charge  amplifier  in  order  to 
prevent  an  overshoot  of  the  signal. 

3.4  Test  Conditions  and  Procedure.  The  four  model  configurations  were  tested  at 
seven  driver  pressure  levels;  namely,  at  P4j  =  P^/P,  -  4,  14,  38,  55,  80,  110,  and 
188.  To  test  the  functionality  of  the  modified  shock  tube  and  to  get  a  feel  for  the 
starting  process  of  the  various  nozzles,  the  models  were  tested  at  the  lowest  speci¬ 
fied  driver  pressure  level  (P41  =  4)  first.  The  tests  then  continued  with  the  next 
highest  driver  pressure  level  (P4J  =  14).  The  operating  driver  pressure  was  1.3  MPa 
(13  bar)  over  atmospheric  pressure,  except  for  the  lowest  driver  pressure  ratio 
where  the  driver  pressure  was  300  kPa  (3  bars)  over  atmospheric  pressure.  The 
actual  initial  test  conditions  are  listed  in  Table  II  for  all  experiments.  The  experi¬ 
mental  driver  pressure  p4  given  in  the  table  is  the  fill  pressure  above  atmospheric 
pressure  and  is  related  to  the  absolute  driver  pressure  by  p4  =  P4  -  P0.  The  indivi¬ 
dual  letters  B,  C,  and  D  indicate  the  window  location  at  which  the  optical  records 
were  taken  (see  Figure  6a).  The  combination  of  B/C  and  C/D  indicates  that  the  win¬ 
dow  was  located  halfway  in  between  the  positions  shown  in  Figure  6a. 

The  test  preparation  was  begun  by  installing  the  desired  model  configuration 
and  the  diaphragm.  The  shock  tube  was  then  closed,  the  Cranz-Schardin  camera 
and  the  pressure  recorders  were  readied,  and  all  electrical  connections  were  checked. 
The  recording  film  holder  was  opened  in  the  darkened  room  and  the  firing  and 
recording  controls  were  activated  before  the  driver  was  pressurized.  Once  the 
driver  had  reached  the  specified  pressure,  the  driven  section  was  evacuated  as  neces¬ 
sary  until  the  desired  driver  pressure  ratio  was  established  (see  Table  II).  The  rup¬ 
turing  of  the  diaphragm  was  initiated  with  an  electric  pulse  of  10-20  kV  to  explode 
the  wire  embedded  between  the  two  Ultraphan  layers  of  the  diaphragm.  The  same 
pulse  triggered  the  camera  and  the  other  recording  equipment  after  some  preset 
delay  time. 
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Config/P41 

TABLE  II.  INITIAL  TEST  CONDITIONS 

Shot#  Window  P0  T0  p4  Pi 

Pos.  mbar  °C  bar  mbar 

P41 

90°  I 

lalf-Nozzle  (N90) 

N90/004 

13249 

B 

990 

22.7 

3.0 

990 

4.030 

it 

13254 

B 

989 

23.5 

it 

989 

•4.033 

N90/014 

13258 

B 

988 

22.6 

13.0 

988 

14.16 

N90/038 

13310 

B 

996 

23.2 

13.0 

368 

38.03 

N90/055 

13301 

B 

992 

24.2 

13.0 

254 

55.09 

it 

13302 

B 

992 

23.7 

ii 

■ 

55.09 

N90/080 

13311 

B 

996 

23.4 

13.0 

174 

80.44 

N90/110 

13314 

B 

993 

23.5 

13.0 

127 

110.2 

N90/188 

13312 

B 

994 

23.9 

13.0 

74 

189.1 

ft 

13313 

B 

993 

23.2 

ii 

it 

189.1 

45°  I 

lalf-Nozzle  (N45) 

N45/004 

13244 

B 

990 

24.6 

3.0 

990 

4.030 

II 

13253 

B 

990 

23.2 

ii 

ti 

4.030 

N45/014 

13259 

B 

988 

23.3 

13.0 

988 

14.16 

H 

13265 

B/C 

995 

24.3 

ii 

995 

14.07 

N45/038 

13309 

B 

995 

22.5 

13.0 

368 

38.03 

N45/055 

13303 

B 

992 

24.3 

13.0 

254 

55.09 

it 

13321 

B/C 

993 

23.8 

it 

ii 

55.09 

N45/080 

13308 

B 

990 

24.2 

13.0 

174 

80.40 

n 

13319 

B/C 

994 

22.3 

ti 

ii 

80.43 

N45/110 

13306 

B 

990 

24.0 

13.0 

127 

110.2 

it 

13320 

B/C 

993 

23.5 

ii 

ii 

110.2 

N45/188 

13307 

B 

990 

24.2 

13.0 

74 

189.1 

n 

13318 

B/C 

990 

24.8 

ii 

ii 

189.1 

The  principle  of  the  electronic  control  circuit  for  the  experimental  process  can 
be  seen  in  Figure  9.  The  preset  delay  time  between  the  ignition  pulse  for  bursting 
the  diaphragm  and  the  first  spark  pulse  of  the  Cranz-Schardin  camera  was  measured 
by  a  time  counter.  Some  interference  from  the  strong  trigger  pulse  was  suffered  on 
the  recording  equipment,  and  a  few  pressure  records  were  lost.  A  test  log  was  kept 
to  record  these  events  together  with  the  initial  test  and  ambient  conditions  for  each 
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shot,  the  trigger  times  for  the  pressure  recorders  and  for  the  Cranz-Schardin  cam 
era,  and  the  spark  frequencies  for  the  camera. 


Config/P41 

rABLE  I 

Shot# 

I.  INITIAL  TEST 

Window  P0 

Pos.  mbar 

CONDITIONS  (continued) 

To  P4  Pi 

°C  bar  mbar 

P4i 

16° 

Half-Nozzle  (N16) 

N 16/004 

13245 

B 

986 

24.7 

3.0 

986 

4.043 

It 

13246 

C 

990 

23.6 

it 

990 

4.030 

N 16/014 

13260 

B 

988 

23.5 

13.0 

988 

14.16 

it 

13264 

C 

994 

23.8 

it 

994 

14.08 

N 16/038 

13317 

B 

990 

24.3 

13.0 

368 

38.02 

it 

13322 

C 

993 

22.2 

ii 

it 

38.02 

N 16/055 

13304 

B 

992 

24.8 

13.0 

254 

55.09 

it 

13323 

C 

993 

22.8 

it 

ll 

55.09 

N 16/080 

13316 

B 

993 

24.0 

13.0 

174 

80.42 

t» 

13324 

C 

993 

23.0 

11 

ll 

80.42 

N 16/1 10 

13305 

B 

990 

23.4 

13.0 

127 

110.2 

it 

13325 

C 

993 

23.7 

it 

it 

110.2 

N 16/188 

13315 

B 

993 

24.1 

13.0 

74 

189.1 

H 

13326 

C 

993 

24.0 

ii 

ll 

189.1 

6° 

Half-Nozzle  (N06) 

N  06/004 

13247 

B/C 

987 

24.0 

3.0 

987 

4.040 

II 

13248 

C/D 

985 

24.5 

ii 

985 

4  046 

N06/014 

13261 

B 

990 

23.6 

13.0 

990 

14.13 

II 

13262 

C/D 

994 

23.0 

ii 

994 

14.08 

it 

13263 

C 

994 

23.4 

" 

994 

14.08 

N06/038 

13327 

C/D 

993 

23.0 

13.0 

368 

38.02 

N 06/05 5 

13328 

C/D 

993 

23.2 

13.0 

254 

55.09 

N06/080 

13329 

C/D 

992 

24.5 

13.0 

174 

80.41 

N06/110 

13330 

C/D 

992 

25.5 

13.0 

127 

110.2 

N06/188 

13331 

C/D 

992 

24.5 

13.0 

74 

189.1 

4.  EVALUATION  OF  EXPERIMENTAL  DATA 


The  nozzle  flow  start-up  phenomena  observed  and  the  data  reduction  methods 
used  are  described  in  this  section.  With  each  test  shot,  24  shadow  or  schlieren  pic¬ 
tures  were  recorded  on  film  by  the  Cranz-Schardin  camera  at  constant  time  inter¬ 
vals.  A  set  of  the  optical  records  was  forwarded  by  EMI  to  BRL  for  data  reduction 
and  analysis.  The  complete  set  of  X-t  diagrams  of  the  evaluated  optical  records  is 
presented  in  Appendix  A. 

The  static  pressure  was  recorded  at  three  available  gage  locations  in  the  tunnel 
ceiling,  270mm,  760mm,  and  889mm  downstream  from  the  center  of  the  test  sec¬ 
tion  window.  The  pressure  was  recorded  at  all  seven  driver  pressure  ratios  for  the 
90°,  45°,  and  16°  divergent  half-nozzles  in  the  window  position  B  for  which  the 
diaphragm  is  located  at  the  upstream  edge  of  the  test  section  window,  100mm  from 
its  center.  A  set  of  the  pressure  records  was  also  forwarded  by  EMI  to  BRL  for 
evaluation. 

4.1  Description  of  Typical  Flow  Phenomena.  The  16°  divergent  half-nozzle 
appears  to  be  best  suited  for  an  explanation  of  the  typical  flow  start-up  phenomena 
observed  in  this  study  because  the  available  field  of  view  (400mm)  is  the  longest 
among  the  four  configurations  under  investigation.  Two  test  runs  "'ere  needed  to 
record  the  flow  throughout  the  entire  length  of  the  16°  half-nozzle  because  the  field 
of  view  of  the  optical  window  was  limited  to  200mm.  One  test  run  covered  the  first 
200  mm  of  the  16°  half-nozzle  (window  position  B),  and  another  test  run  covered  the 
range  from  200 -400 mm  behind  the  diaphragm  location  (window  position  C). 
Therefore,  each  snapshot  in  time  as,  e.g.,  shown  in  Figure  10  is  composed  of  two 
separate  photographic  records.  A  heavy  vertical  line  halfway  through  the  nozzle  of 
some  of  the  snapshots  indicates  that  the  fields  of  view  of  the  two  photographs  did 
not  quite  overlap. 

A  low-pressure  case  was  chosen  to  illustrate  the  flow  start-up  phenomena 
because  all  the  phenomena  of  interest  remain  within  the  field  of  view  until  the  start¬ 
up  process  of  the  nozzle  is  completed.  A  sample  sequence  of  shadowgraphs  for  the 
16°  divergent  half-nozzle  (N16)  at  the  driver  pressure  ratio  P4i  =  14  (Code  N 16/014) 
is  presented  in  Figures  10a  and  10c.  The  major  nozzle  flow  phenomena  were  traced 
from  these  photographs  and  are  illustrated  in  Figures  10b  and  lOd.  Figures  10a 
and  10b  show  the  early  phase  of  the  flow  start-up  process,  and  Figures  10c  and  lOd 
show  the  late  phase  cf  the  flow  start-up  process.  Flow  phenomena  like  the  incident 
shock,  the  contact  surface,  boundary  layer  build-up,  regions  of  flow  separation,  and 
a  system  of  multiple  recompression  shocks  can  be  seen  developing  and  moving 
within  and  past  the  field  of  view.  They  are  described  in  the  following  paragraphs. 

The  incident  shock  begins  to  form  immediately  after  the  flow  is  initiated  by 
bursting  the  diaphragm.  Since  the  diaphragm  breaks  into  fragments  of  uneven  sizes, 
the  opening  process  is  asymmetrical  and  causes  an  uneven  flow  progress  across  the 
height  and  depth  of  the  throat.  Initially,  the  shock  consists  of  a  system  of  two  to 
three  curved  compression  waves  which  catch  up  with  each  other  inside  the  throat 
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channel  (Figures  10a, b  at  106  (is  after  diaphragm  break)  and  reflect  from  the  ceiling 
and  the  floor  of  the  throat  channel  combining  eventually  into  one  incident  shock 
(206 (is).  This  formation  period  may  take  up  to  about  500  (is  at  the  low  end  of  the 
driver  pressure  range,  by  which  time  the  shock  has  progressed  well  into  the  diver¬ 
gent  section  of  the  half  nozzle  (606 (is).  A  shock  reflection  is  generated  at  the  ceil¬ 
ing  of  the  half  nozzle  exit  (Figures  10c, d  at  821  (is)  proceeding  across  the  flow  area 
and  reflecting  later  on  at  the  floor  plate  of  the  expansion  channel. 

The  contact  surface  between  the  expanding  driver  gas  and  the  gas  in  the  low- 
pressure  region  of  the  divergent  half-nozzle  follows  the  incident  shock  at  close 
range.  It  can  be  easily  detected  because  the  associated  flow  turbulence,  which  it 
brings  with  it,  changes  the  texture  of  the  background  in  the  photographs  (see 
Figures  10a  and  10c).  Because  of  the  turbulence,  the  contact  front  is  very  rugged. 
Often,  darker  strings  perpendicular  to  the  flow  direction  can  be  detected  in  the  tur¬ 
bulent  region  behind  the  contact  front  (e.g.,  Figure  10a,  606(is).  It  is  thought  that 
these  strings  are  combustion  products  of  the  exploding  wire  that  ruptures  the 
diaphragm,  rather  than  compression  waves  which  dissipate  with  time  (Figure  10c, 
1.321ms).  Once  the  contact  surface  has  entered  the  divergent  section  of  the  half¬ 
nozzle,  a  variety  of  phenomena  begin  to  occur  simultaneously. 

An  expansion  fan  begins  to  develop  at  the  exit  of  the  throat  channel,  at  the 
corner  with  the  divergent  section  of  the  half  nozzle.  It  is  best  seen  in  Figure  10c;  in 
Figure  10a,  the  diaphragm  fragments  exiting  the  throat  channel  obscure  the  image  of 
the  developing  expansion  fan.  Inside  the  divergent  section  of  the  half-nozzle,  a 
recompression  shock  begins  to  form  (606 (is).  This  shock  develops  into  a  bottle 
shock  spawning  a  second  shock  ahead  of  itself  (Figure  10c,  1.322ms).  The  second 
shock  also  develops  into  a  bottle  shock  and  may  spawn  a  third  shock  ahead  of  itself 
in  the  process  (1.822  ms).  Extensive  flow  separation  regions  at  both  ceiling  and 
floor  of  the  divergent  half-nozzle  develop  during  this  process  due  to  the  boundary- 
layer  build-up  in  the  flow  channel.  With  the  boundary  layer  well  developed,  the 
system  of  bottle  shocks  changes  into  a  system  of  multiple  X-type  recompression 
shocks.  The  start-up  process  of  the  nozzle  comes  to  completion  when  this  system  of 
multiple  X-type  shocks  becomes  stationary  either  inside  the  nozzle,  as  in  the  present 
case  shown  in  Figures  10c  and  lOd  (2.322ms),  or  downstream  from  the  divergent 
half-nozzle  for  high  driver  pressure  ratios,  as  suggested  in  the  photographs  shown  in 
Figures  lib  and  12. 

Throat  shocks  appear  inside  the  throat  channel  as  a  consequence  of  the 
boundary-layer  build-up  (406 (is  and  later).  As  many  as  six  shocks  have  been 
observed  at  the  lowest  pressure  level  (P41  =  4).  Two  of  these  warrant  special  men¬ 
tion  because  of  their  general  presence  at  almost  all  test  conditions.  The  narrowest 
cross  section  exists  at  the  diaphragm  location.  Once  the  flow  reaches  the  sonic  con¬ 
dition  in  this  cross  section,  it  expands  supersonically  upon  exiting  into  the  throat 
channel  and  a  standing  X-type  recompression  shock  appears  between  10-15  mm 
downstream  from  the  diaphragm  location  (606 (is)  at  all  pressure  levels  ex'ept  the 
lowest  one.  A  second,  bifurcated  recompression  shock  appears  in  the  throat  channel 
at  about  50-60 mm  downstream  from  the  diaphragm  location  and  remains  there 
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throughout  the  start-up  process  of  the  nozzle  (821  (is  and  later).  A  strong  shock  at 
first,  it  separates  into  multiple  bifurcated  or  multiple  X-type  shocks  at  later  times. 

Perhaps  the  biggest  surprise  in  these  photographs  is  the  discovery  of  the  asym¬ 
metric  flow  through  the  divergent  half-nozzle.  It  was  the  expectation  of  the  experi¬ 
menters  to  obtain  one-half  of  a  symmetrical  flow  pattern  above  the  center  axis  of  the 
nozzle.  This  was  accomplished  for  the  flow  ahead  of  and  up  to  the  contact  surface. 
Behind  the  contact  surface,  however,  the  boundary-layer  build-up  along  the  floor 
plate  and  the  ceiling  eventually  causes  the  flow  to  separate.  The  flow  then  seeks  its 
own  path  by  turning  and  more  or  less  following  the  centerline  of  the  half-nozzle. 

4.2  Vortices.  Fragments,  and  Other  Curiosities.  The  90°  half-nozzle  is  a  nozzle 
without  a  divergent  section.  Two  examples  of  the  flow  start-up  process  in  the  90° 
half-nozzle  are  presented  in  Figure  11.  Figure  11a  shows  a  sequence  of  selected  sha¬ 
dowgraphs  for  the  lowest  driver  pressure  ratio  tested,  P4j  =  4  (Code  N90/004). 
Recompression  shocks  can  be  seen  in  the  throat  and  in  the  flow  channel  downstream 
from  the  throat  exit  where  the  flow  begins  to  turn,  but  the  typical  X-type  shock  near 
the  diaphragm  is  not  seen  at  this  pressure  level.  A  system  of  five  recompression 
shocks  develops  in  the  throat  of  which  one  shock  becomes  stationary  at  the  throat 
exit  and  another  one  becomes  stationary  a  short  distance  downstream  from  the 
throat  exit. 

The  interesting  phenomenon  of  this  case  (N90/004)  is  the  big  vortex  that  fills 
the  corner  downstream  from  and  above  the  exit  from  the  throat  channel  while  the 
flow  continues  straight  ahead  in  a  channel  the  height  of  the  throat  along  the  floor  of 
the  expansion  tunnel.  The  vortex  can  be  seen  beginning  to  develop  in  Figure  11a  at 
0.302ms  and  filling  the  entire  dead-water  region  at  about  1.302ms.  This  vortex 
flow  suggests  that  a  sizable  amount  of  internal  energy  is  used  up  by  the  filling  of  the 
dead-water  region  formed  by  the  area  discontinuity  at  the  throat  exit;  and  this 
energy  loss  should  have  an  influence  on  the  downstream  pressure  recovery. 

This  pressure  level  is  one  of  two  at  which  such  a  big  vortex  was  observed  for 
the  90°  half-nozzle  configuration.  At  the  next  higher  pressure  level,  P4i  =  14,  the 
vortex  is  equally  big  but  fills  the  dead-water  region  faster  and  takes  on  a  different 
appearance  because  the  structure  of  the  turning  flow  is  distinctly  different  from  the 
present  case.  The  lack  of  observing  such  distinct  vortices  at  higher  pressure  levels 
(i.e.,  P4j  =  38,  and  above)  does  not,  however,  imply  that  they  do  not  exist.  Rather 
it  is  surmised  that  due  to  the  turning  of  the  flow  and  the  faster  equalization  of  the 
static  pressure  in  the  corner,  they  may  have  been  missed  at  the  relatively  coarse 
sampling  sequence  employed  in  these  tests. 

The  mainstream  flow  out  of  the  throat  in  Figure  11a  continues  straight  ahead 
for  the  entire  length  of  the  test  section  without  turning  (see  picture  at  2,102ms)  and 
the  schlieren  picture  at  1.702ms  shows  a  very  distinct  slip  stream  commencing  at  the 
throat  exit  and  continuing  through  the  range  of  observation.  Unfortunately,  the 
flow  further  downstream  from  the  throat  exit  could  not  be  observed  because  it  was 
not  possible  to  mount  this  configuration  further  upstream  from  the  test  section 
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window.  This  proved  a  major  handicap  at  the  higher  pressure  ratios  where  the  flow 
phenomena  of  interest  all  moved  past  the  field  of  view  before  the  start-up  process 
was  completed. 

Figure  11b  shows  a  sequence  of  selected  shadowgraphs  for  the  start-up  process 
at  the  highest  driver  pressure  ratio  investigated,  P41  =  188  (Code  N90/188).  A  bottle 
shock  can  be  seen  to  develop  early  in  the  start-up  process  (240-390p.s)  before  the 
recompression  shock  moves  out  of  the  field  of  view.  From  240 p.s  on,  the  flow 
appears  to  expand  linearly  at  an  increasing  angle  and  steadying  at  45°  after  690 jxs. 
The  flow  expansion  appears  to  be  disturbed  by  the  strong  reflection  of  the  incident 
shock  from  the  ceiling  at  540  ns.  By  this  time,  the  recompression  shock  has  moved 
beyond  ihe  field  of  view,  and  the  start-up  process  has  escaped  observation  to  its 
completion. 

The  flow  start-up  process  in  the  45°  half-nozzle  is  illustrated  in  Figure  12  for 
the  driver  pressure  ratio,  P4i  =  80  (Code  N45/080).  For  this  configuration,  a  field  of 
view  of  300mm  was  available  which  allowed  us  to  observe  the  start-up  How 
phenomena  for  a  longer  distance  downstream  from  the  throat  exit  than  for  the  90° 
half-nozzle.  Still,  the  field  of  view  is  too  limited  to  record  all  recompression  shocks 
up  to  the  completion  of  the  start-up  period.  The  shadowgraph  sequence  in  Figure 
12a  shows  that  a  bottle  shock  has  formed  at  385 p.s,  before  the  diaphragm  fragments 
enter  the  divergent  part  of  the  half-nozzle.  At  later  times  (Figure  12b),  the  expand¬ 
ing  flow  and  the  bottle  shock  are  disturbed  by  the  reflected  shock  from  the  ceiling 
of  the  expansion  tunnel  and  by  the  diaphragm  fragments. 

The  fragments  of  the  ruptured  diaphragms  present  an  unavoidable  disruption 
of  the  flow  field  in  all  test  cases.  In  the  example  shown  (Figure  12b),  the 
diaphragm  fragments  are  seen  to  disrupt  the  bottle  shock  from  803  p.s  to  1.103  ms. 
As  the  fragments  leave  the  scene,  the  bottle  shock  recovers  and  become*  stationary 
within  the  field  of  view  (1.463  ms).  However,  the  second  recompression  shock 
which  had  already  formed  ahead  of  the  bottle  shock  at  803 p.s  has  moved  beyond  it. 
The  throat  shocks  also  are  clearly  visible  in  all  the  photographs  of  Figure  12. 

The  start-up  process  of  the  nozzle  is  also  characterized  by  the  flow  separation 
along  the  upper  wall  of  the  divergent  nozzle.  Initially,  the  flow  separates  at  the 
throat  exit.  A  distinct  slip  stream  can  be  seen  at  385 jxs  commencing  at  the  throat 
exit.  As  the  start-up  process  continues,  the  separation  point  moves  up  along  the 
wall  toward  the  nozzle  exit  at  the  ceiling  of  the  expansion  channel  until  a  stable  con¬ 
dition  is  reached  (1.103ms).  The  dead-water  region  near  the  corner  is  considerably 
smaller  than  the  one  observed  with  the  90°  half-nozzle  and  much  of  the  flow  vorti- 
city  is  avoided  with  the  addition  of  the  45°  divergent  nozzle.  The  significance  of  the 
effect  which  the  45°  divergence  has  on  the  flow  will  be  discussed  in  the  following 
section. 

4.3  Data  Reduction  Procedures.  The  shadow  and  schlieren  pictures  were 
evaluated  by  measuring  the  distances  of  the  various  flow  phenomena  from  a  refer¬ 
ence  point  in  the  photograph  with  the  help  of  a  digitizing  pad  that  was  controlled  by 
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a  Zenith  PC.  The  measured  locations  of  these  flow  phenomena,  e.g.,  incident 
shock,  contact  surface,  various  recompression  shocks,  and  some  diaphragm  frag¬ 
ments,  with  reference  to  the  diaphragm  location  (X  =  0),  were  then  recorded  as 
functions  of  time  measured  from  the  bursting  of  the  diaphragm  (t  =  0).  From  these 
records,  X-t  diagrams  were  generated  on  an  Alliant  FX/1  computer  using  PVI  (Pre¬ 
cision  Visuals,  Inc.  1986)  software.  The  data  points  for  the  window  positions,  B/C 
(100-300mm),  C  (200-400mm),  and  C/D  (300-500mm),  downstream  from  the 
diaphragm  location  were  matched  with  the  data  points  for  the  window  position  B 
(0-200 mm)  by  adding  a  time  difference  to  the  24  recorded  spark  times  so  that  the 
data  points  for  the  incident  shock  recorded  in  the  downstream  window  positions  C 
and  D  would  continue  on  a  straight  line  drawn  through  the  last  three  data  points  for 
the  incident  shock  recorded  in  window  position  B.  All  X-t  diagrams  are  shown  in 
Appendix  A. 

It  was  learned  from  the  X-t  diagrams  that  the  data  points  for  the  incident  shock 
generally  follow  a  cubic  polynomial.  From  the  X-t  diagram  for  case  N 16/014,  e.g., 
one  can  learn  that  the  incident  shock  and  the  contact  surface  initially  accelerate  in 
the  throat  of  the  nozzle,  become  constant  for  a  short  period  of  time,  and,  finally, 
decelerate  in  the  divergent  section  of  the  nozzle.  Therefore,  a  least-sqares  polyno¬ 
mial  of  third  degree  was  fitted  to  the  data  for  the  incident  shock  to  determine  the 
intercept  in  time  and  this  value  was  then  used  to  correct  the  spark  times  so  that  the 
data  for  the  incident  shock  and  the  contact  surface  would  pass  through  the  origin  of 
the  plot. 

Measurements  of  the  static  pressure  were  taken  at  three  available  locations 
along  the  ceiling  of  the  expansion  tunnel,  370mm,  860mm,  and  989mm  down¬ 
stream  from  the  diaphragm  location.  Examples  are  shown  in  Figures  13-15.  Noise 
caused  by  the  high-voltage  trigger  signal  can  be  detected  at  time  zero.  Most  of  the 
pressure  records  show  two,  and  sometimes  three,  shocks  arriving  at  the  gage  posi¬ 
tions  in  short  succession  of  each  other  (about  10-50ps  apart).  In  a  first  attempt  at 
evaluation,  these  pressure  peaks  were  measured,  scaled,  and  recorded  with  their 
respective  arrival  times  for  each  of  the  three  gages.  A  comparison  of  the  shock 
recovery  for  the  three  half-nozzles  using  the  pressure  peaks  of  the  pressure  his¬ 
tories,  however,  proved  unsatisfactory  and  was  discarded.  Consequently,  it  was 
decided  to  use  a  different  approach. 

In  a  second  attempt  at  evaluation,  the  height  of  the  pressure  plateau  behind  the 
incident  shock(s)  was  determined  from  the  pressure  records  averaging  out  the  pres¬ 
sure  fluctuations  by  inspection.  This  judgement  appeared  somewhat  subjective 
because  the  fluctuations  of  the  pressure  traces  are  considerable  indeed;  however,  this 
approach  was  more  successful  than  the  first  one.  In  one-third  of  the  records,  the 
pressure  plateau  agreed  with  the  first  pressure  peak  (Figure  13),  in  another  third  it 
coincided  with  the  second  pressure  peak  (Figure  14),  and  in  the  last  third  of  the 
records  it  fell  in  between  the  two  pressure  peaks  (Figure  15).  From  the  measure¬ 
ments,  the  shock  strength,  P2i,  was  determined  and  plotted  versus  the  driver  pres¬ 
sure  ratio,  P4i,  and  the  shock  Mach  number,  M,.  The  results  are  discussed  in  the 
following  section. 
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5.  DISCUSSION  OF  RESULTS 


The  objective  of  this  experimental  investigation  was  to  gather  information  on 
the  flow  start-up  phenomena  in  convergent-divergent  nozzles  such  as  are  used  in 
large  blast  simulators.  Of  particular  interest  was  the  time  required  for  the  incident 
shock  to  form,  the  influence  of  the  divergent  nozzle  angle  on  the  shock  recovery, 
the  flow  start-up  period,  and  a  comparison  with  computational  predictions  of  the 
flow  start-up  process.  Subsequently,  the  shadowgraphs,  X-t  diagrams,  and  pressure 
records  were  evaluated  with  this  objective  in  mind.  The  flow  start-up  period,  the 
shock  formation  process,  the  flow  expansion  in  the  divergent  nozzles,  and  the  pres¬ 
sure  recovery  for  the  16°,  45°,  and  90°  nozzles  were  summarized  in  graphs  of  the 
defined  nozzle  characteristics  versus  driver  pressure  ratio  or  the  shock  Mach  number 
and  are  discussed  in  this  section. 

The  6°  half-nozzle  was  excluded  from  the  discussion  because  the  available  data 
were  too  limited  in  scope.  Except  at  driver  pressure  ratios  P4i  =  4  and  14,  optical 
records  were  obtained  only  in  window  position  C/D  which  afforded  a  view  of  the 
area  discontinuity  at  the  nozzle  exit.  The  area  discontinuity,  it  turned  out,  had  an 
effect  comparable  to  that  of  the  90°  nozzle  in  that  the  nozzle  acted  like  an  extended 
throat  channel  and,  with  the  exception  of  the  formation  of  the  incident  shock  and 
the  throat  shocks,  the  start-up  phenomena  of  interest  occured  mostly  downstream 
from  the  nozzle  exit.  The  field  of  view  in  this  position  proved  to  be  too  limited  to 
successfully  record  the  flow  start-up  processes. 

5.1  Flow  Start-up  Period.  The  start-up  process  of  the  nozzle  was  defined  in  Sec¬ 
tion  2.3  as  the  time  from  the  rupturing  of  the  diaphragm  to  the  establishment  of  a 
quasisteady  flow  pattern  in  the  nozzle.  All  flow  phenomena  observed  during  the 
start-up  process  of  the  half-nozzles  are  well  known  from  previous  research  of  the 
steady  flow  in  ducts  and  nozzles  at  various  exit  pressure  ratios  of  the  nozzle 
(Emmons  1958).  The  only  difference  for  unsteady  flow  appears  to  be  that  all  these 
various  phenomena  appear  together  in  the  course  of  one  flow  sequence  changing 
from  one  into  another  with  time.  For  example,  a  single,  strong,  bifurcated  throat 
shock  changes  with  time  into  weaker,  multiple  bifurcated  shocks  and  still  later  into 
multiple  X-type  shocks  in  response  to  the  changing  back  pressure  ratio  of  the  noz¬ 
zle.  This  holds  also  true  for  the  system  of  recompression  shocks  which  develops  in 
the  divergent  nozzle  and  downstream  from  its  exit,  respectively. 

The  flow  start-up  period  was  determined  from  the  X-t  diagrams  for  three  half¬ 
nozzles.  It  is  shown  in  Figure  16  as  a  function  of  the  driver  pressure  ratio.  The 
data  show  that  the  flow  start-up  period  decreases  for  all  nozzles  with  increasing 
divergent  angle.  The  flow  start-up  period  for  the  16°  half-nozzle  could  be  deter¬ 
mined  for  three  pressure  levels  only  (P41  =  14,  38,  55)  because  either  the  recording¬ 
time  periods,  or  the  field  of  view  for  the  other  test  cases  were  insufficient.  These 
data  points  are  close  above  those  for  the  45°  nozzle.  The  relatively  large  difference 
in  start-up  time  for  the  SO0  half-nozzle  as  compared  to  the  45°  half-nozzle  can  be 
explained  by  the  fact  that  the  flow  start-up  period  for  this  nozzle  was  determined 
from  the  data  for  the  throat  shocks  because  the  field  of  view  behind  the  throat  exit 
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was  too  limited  to  observe  the  entire  start-up  process.  It  is  surmised  that  if  the 
entire  start-up  process  could  have  been  observed,  the  differences  in  start-up  time 
between  the  45°  and  90°  half-nozzles  would  have  been  found  smaller,  because  in 
those  cases  in  which  it  could  be  observed  (e.g.,  case  N16/014  and  N45/080),  the  X-t 
diagrams  reveal  that  the  system  of  bottle  shocks  downstream  from  the  nozzle  exit 
becomes  stationary  at  somewhat  later  times  than  the  throat  shocks. 

5.2  Shock  Formation.  In  the  general  description  of  the  flow  phenomena  in  Sec¬ 
tion  4.1,  it  was  mentioned  that  the  incident  shock  starts  out  as  a  series  of  compres¬ 
sion  waves  which  catch  up  with  each  other  in  the  throat  channel  of  the  nozzle  (e.g., 
Figure  12,  75  (is).  The  incident  shock  is  further  built-up  by  shock  reflections  from 
the  bottom  or  the  ceiling  of  the  throat  channel  (Figure  12,  135 (is)  and  appears  to  be 
well-formed  a  short  time  after  it  has  emerged  from  the  throat  channel  (Figure  12, 
255 (is).  These  observations  could  be  made  in  all  test  cases  and  Figures  10-12 
present  typical  examples  of  the  formation  of  the  incident  shock.  The  approximate 
period  and  distance  of  shock  formation  were  determined  by  inspection  from  the  opt¬ 
ical  records  for  each  test  case  and  the  results  are  plotted  in  Figure  17. 

Figure  17a  shows  the  time  it  takes  the  incident  shock  to  be  well-formed  as  a 
function  of  the  driver  pressure  ratio  for  the  three  half-nozzles.  The  results  indicate 
that  the  incident  shock  forms  fastest  in  the  16°  half-nozzle  and  slowest  in  the  90° 
half-nozzle.  The  time  difference  between  the  two  configurarions  appears  to  be 
approximately  50 (is  over  the  whole  range  of  driver  pressure  ratios.  The  shock  for¬ 
mation  in  the  45°  half-nozzle  presents  a  transition  between  the  other  two  nozzles 
starting  out  with  periods  even  longer  than  in  the  90°  nozzle  at  the  two  lowest  pres¬ 
sure  levels,  and  ending  up  with  shock  formation  periods  almost  identical  to  the  16° 
nozzle  at  pressure  ratios  above  40.  Figure  17b  presents  the  interesting  result  that 
the  shock  formation  distance  is  practically  constant  over  the  whole  pressure  range 
tested. 

The  differences  found  in  the  shock  formation  periods  for  the  three  nozzle 
angles  may  be  significant  for  the  LB/TS  design  if  these  times  scale  with  the  volume 
of  the  facility.  The  statistical  significance  of  the  differences  found  is  minimized 
however  by  the  tolerances  involved  in  their  determination.  The  tolerance  associated 
with  each  data  value  is  proportionate  to  the  time  interval  which  was  used  for 
sequencing  the  light  sparks  of  the  Cranz-Schardin  camera.  At  the  lower  end  of  the 
pressure  range,  the  time  interval  between  data  points  is  100 (is,  and  at  the  upper  end 
of  the  pressure  range  it  narrows  down  to  40 (is.  As  a  result,  the  differences  in 
shock  formation  period  between  the  three  nozzles  may  be  insignificant  altogether 
because  they  are  of  the  same  magnitude  as  the  tolerances  involved  in  measuring 
them. 

5.3  Flow  Expansion  and  Separation.  The  most  significant  flow  feature  of  the  90° 
half-nozzle  is  that  the  flow  forms  its  own  divergent  nozzle.  As  the  driver  gas 
behind  the  contact  surface  leaves  the  throat  channel,  the  flow  begins  to  expand  at 
some  angle  depending  on  the  back  pressure  at  the  throat  channel  exit.  This  angle 
can  be  identified  by  the  slip  stream  visible  in  the  shadowgraphs  (e.g.,  Figure  11). 
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The  expansion  angle  approaches  a  quasi-steady  value  once  the  local  disturbances 
caused  by  the  shock  reflections  from  the  ceiling  and  the  floor  have  died  down  and 
the  diaphragm  fragments  have  passed  (e.g.,  Figure  lib,  915  p.s).  The  flow  expan¬ 
sion  angle  was  determined  by  reading  and  averaging  the  maximum  and  minimum 
angles  of  the  slip  stream  in  the  shadowgraphs  and  schlieren  pictures  at  late  times 
close  to  the  completion  of  the  start-up  process.  It  is  plotted  versus  the  driver  pres¬ 
sure  ratio  in  Figure  18a  and  compared  to  the  flow  expansion  angles  in  the  45°  and 
16°  half-nozzles  which  were  determined  in  the  same  manner.  The  results  in 
Figure  18a  show  that  the  flow  stays  attached  in  the  45°  nozzle  for  driver  pressure 
ratios  of  55  and  above  but  seeks  its  own  expansion  angle  at  lower  pressure  ratios. 
In  the  16°  half-nozzle,  the  flow  stays  attached  for  driver  pressure  ratios  above  14. 

The  complimentary  phenomenon  to  flow  expansion  is  the  separation  of  the 
flow  from  the  walls  of  the  diverging  duct  that  is  formed  by  the  nozzle  configura¬ 
tions.  The  flow  separation  angle  for  the  three  half-nozzles  is  plotted  versus  the 
driver  pressure  ratio  in  Figure  18b.  In  the  90°  half-nozzle,  the  flow  has  no  choice 
but  to  separate  upon  exiting  at  the  exit  of  the  throat  channel.  The  flow  separation 
angle  is  defined  in  the  same  manner  as  the  flow  expansion  angle  and  is  therefore 
equal  to  the  latter.  In  the  16°  nozzle,  on  the  other  hand,  flow  separation  is  almost 
nonexistent.  Only  at  the  two  lowest  driver  pressure  ratios  does  the  flow  separate 
halfway  through  the  nozzle  due  to  the  boundary  layer  build-up.  At  pressure  levels 
above  14,  the  flow  stays  attached  to  the  nozzle  walls  almost  to  the  nozzle  exit,  where 
it  forms  a  small  dead-water  region  to  negotiate  the  corner  at  the  ceiling  of  the 
expansion  tube. 

In  the  45°  divergent  half-nozzle,  flow  separation  occurs  at  the  three  lowest 
driver  pressure  levels  tested,  i.e.,  P41  =  4,  14,  and  38.  At  these  pressures,  the  flow 
separates  at  the  exit  corner  of  the  throat  channel  for  the  entire  flow  start-up  period. 
At  higher  driver  pressure  levels,  flow  separation  occurs  early  on  in  the  start-up  pro¬ 
cess  at  the  exit  corner  of  the  throat  channel  but  as  the  start-up  process  continues,  the 
separation  point  moves  downstream  along  the  ceiling  toward  the  exit  of  the  diver¬ 
gent  nozzle.  Toward  the  end  of  the  start-up  process,  the  separation  point  becomes 
stationary  at  some  point  along  the  ceiling  of  the  divergent  nozzle.  The  flow  then 
stays  attached  to  the  ceiling  of  the  divergent  section  of  the  half-nozzle. 

The  flow  separates  at  a  point  near  the  corner  of  the  nozzle  exit,  creating  a 
dead-water  region  to  negotiate  the  corner  formed  by  the  divergent  nozzle  and  the 
expansion  tube.  The  function  of  the  flow  separation  angle  for  the  45°  half  nozzle  in 
Figure  18b  shows  a  maximum  about  the  driver  pressure  ratio,  P41  =  40,  and  slightly 
smaller  separation  angles  for  the  higher  pressure  ratios.  It  is  noted  that  this  is  not 
the  same  angle  as  the  flew  expansion  angle  shown  in  Figure  18a.  In  comparison  to 
the  90°  half-nozzle,  however,  the  flow  through  the  45°  half-nozzle  showing  separa¬ 
tion  angles  of  the  same  magnitude  avoids  much  of  the  vorticity  found  in  the  flow 
field  at  the  area  discontinuity  of  the  90°  half-nozzle. 
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5.4  Pressure  Recovery.  The  static  pressure  history  was  recorded  by  pressure 
gages  mounted  in  three  locations  along  the  ceiling  of  the  expansion  tunnel,  i.e., 
370mm  (gage  #37),  860mm  (gage  #  86),  and  989mm  (gage  #  99)  downstream 
from  the  diaphragm  location.  The  incremental  shocks  seen  in  these  records  follow¬ 
ing  the  incident  shock  are  most  likely  caused  by  reflections  of  the  curved  shock  front 
from  the  walls  of  the  expansion  tunnel  as  shown  by  Amann  (1985).  From  the 
records  of  the  pressure  histories  at  these  gage  locations,  the  plateau  pressure  behind 
the  incident  shock  was  determined  and  plotted  versus  (a)  the  driver  pressure  ratio 
and  (b)  the  shock  Mach  number  for  all  half-nozzles  tested.  An  approximate  smooth 
curve  or  a  linear  least-squares  regression  was  fitted  through  the  data  points  for  each 
half-nozzle.  The  results  are  shown  in  Figures  19-21  for  the  three  gage  locations. 

The  results  of  shock  strength,  P2i,  versus  driver  pressure  ratio  P41,  at  gage  #  37 
(Figure  19a),  closest  to  the  diaphragm  location,  do  not  present  conclusive  evidence 
for  a  loss  of  shock  overpressure  as  a  function  of  the  angle  of  divergence  in  the  noz¬ 
zles.  Over  the  whole  range  of  driver  pressure  ratios,  both  the  90°  as  well  as  the  16° 
half-nozzles  show  nearly  identical  shock  overpressures  which  are  higher  than  the 
corresponding  values  for  the  45°  half-nozzle.  However,  when  plotted  versus  the 
shock  Mach  number,  Ms,  (Figure  19b),  the  shock  strength  in  the  16°  half-nozzle 
appears  to  be  higher  than  in  the  other  two  nozzles  for  the  same  shock  Mach 
number.  Therefore,  the  results  are  apparently  somewhat  ambiguous  at  this  gage 
location  because  of  the  closeness  to  the  nozzle  exits.  But,  at  distances  further  away 
from  the  nozzle  exits,  at  gages  #  86  and  #  99,  the  loss  of  shock  overpressure  with 
increasing  angle  of  divergence  is  defined  clearly. 

The  approximate  curve  fits  of  the  shock  strength  versus  driver  pressure  at 
gage  #  86  in  Figure  20a  suggest  that  the  16°  divergent  half-nozzle  produces  an 
enhanced  shock  strength  compared  t .  that  obtained  without  a  divergence,  i.e.,  the 
90°  half-nozzle.  The  difference  in  shock  strength  between  the  two  nozzles  amounts 
to  7.8  percent  at  the  highest  driver  pressure  ratio  (P4i  =  188).  The  gain  in  shock 
strength  behind  the  45°  nozzle  is  3.5  percent  compared  to  the  shock  strength  behind 
the  90°  half-nozzle  at  the  maximum  pressure  ratio.  The  corresponding  data  at 
gage  #99  seem  to  indicate  an  enhancement  of  11  percent  of  the  shock  strength  for 
the  16°  half-nozzle  and  of  7.8  percent  for  the  45°  half-nozzle  over  the  shock  strength 
for  the  90°  half-nozzle  at  the  upper  end  of  the  pressure  scale.  These  results  agree 
with  those  of  small-scale  LBS-model  experiments  carried  out  at  the  BRL  (Coulter 
1987;  Gion,  Private  Communication)  indicating  a  10  percent  enhancement  in  the 
shock  strength  for  drivers  equipped  with  6°  divergent  nozzles. 

Plotting  the  shock  pressure  ratios  against  the  shock  Mach  number 
(Figures  19b -21b),  and  fitting  them  with  a  linear  least-squares  regression,  reveals 
the  uncertainties  involved  in  the  measurements.  Ideally,  all  data  points  would  fall 
on  a  straight  line,  but  in  reality  there  exisf  inaccuracies  in  measurement.  The  corre¬ 
lation  is  very  good,  however;  all  regressions  have  correlation  factors  above  0.9974 
and  the  correlations  improve  with  distance  from  the  diaphragm  location  to  better 
than  0.9992  at  the  gage  location  furthest  away  from  the  diaphragm. 
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An  estimate  of  the  uncertainties  of  the  results  was  also  attempted  by  comparing 
the  data  points  from  gage  #  86  with  those  from  gage  #  99.  The  records  of  the  shock 
overpressure  of  these  gage  positions  indicate  that  the  pressure  loss  between  them  is 
negligible  if  existent  at  all.  The  distance  between  the  two  gages  is  only  129mm, 
and,  in  many  instances,  a  higher  pressure  was  recorded  at  the  more  distant  gage. 
The  relative  differences  in  the  measurements  of  the  shock  pressure  at  these  gages, 
5P86-  99  =  (P86-  P99 )/P99>  may  therefore  be  used  to  determine  the  experimental 
uncertainties.  The  uncertainties  of  the  experimental  measurements  are  defined  by 
the  arithmetic  mean,  x,  of  the  differences  of  all  these  pressure  measurements  and 
their  standard  deviation,  o„_  1.  The  arithmetic  mean  was  determined  to  x  =  0.0175, 
or  1.75  percent  with  a  standard  deviation  of  o„_  1  =  0.0136,  or  1.36  percent  of  the 
relative  differences  in  the  pressure  measurements.  This  result  implies  that  the  loss 
in  shock  strength  between  the  16°  and  the  90°  half-nozzles  of  about  10  percent  is  sta¬ 
tistically  significant  and  therefore  real. 

5.5  Comparison  of  Experimental  Data  with  Computational  Results.  A  systematic 
comparison  of  the  experimental  data  of  this  study  with  BRL-Q1D  code  computa¬ 
tional  results  was  carried  out  and  is  presented  in  Appendix  B.  The  BRL-Q1D  code 
(Opalka  and  Mark  1986)  is  a  computer  program  written  in  FORTRAN  77  language 
that  uses  quasi-one-dimensional,  adiabatic,  inviscid,  numerical  algorithms  to  solve 
the  Euler  equations.  The  code  incorporates  two  computational  techniques.  One  is 
an  implicit  finite-difference  technique  developed  by  Beam  and  Warming,  the  other  is 
an  explicit  finite-difference  technique  according  to  MacC'ormack.  These  techniques 
are  applied  to  the  quasi-one-dimensional  Euler  equations  in  their  weak  conservation 
form.  The  governing  equations  are  derived  for  arbitrary  geometries  and 
transformed  into  a  uniform  computational  grid  retaining  the  weak  conservation 
form.  Linear  and  non-linear  gridding  options  are  available.  The  implicit  technique 
(Beam  and  Warming  1979)  was  used  in  this  application.  Central  spatial  differencing 
casts  the  difference  equations  into  a  block-tridiagonal  structure  which  is  solved  for 
the  increments  in  the  dependent  variables  at  each  successive  time  step. 

The  BRL-Q1D  code  can  provide  only  a  rough,  qualitative  estimate  of  the  com¬ 
plex  flow  start-up  process  in  convergent-divergent  nozzles  because  it  is  limited  to 
one  flow  dimension  and  does  not  model  viscous  effects.  Generally,  the  computa¬ 
tional  results  are  fair  considering  the  inherent  limitations  of  the  code.  The  flow 
phenomena  observed  in  the  experiment  are  generally  reproduced  by  the  computa¬ 
tions  as  long  as  they  are  not  strongly  dependent  upon  viscous  effects.  With  the 
exception  of  the  90°  half-nozzle,  the  experimental  incident  shock  and  contact  surface 
are  matched  by  the  computational  data  for  low  driver  pressures.  With  increasing 
driver  pressure,  however,  the  experimentally  observed  initial  flow  acceleration  and 
shock  formation  phase  becomes  more  pronounced.  Hence,  an  increasing 
discrepancy  between  the  experimental  and  the  computational  paths  for  the  incident 
shock  and  the  contact  surface  is  observed  with  increasing  driver  pressure. 

The  system  of  recompression  shocks  observed  in  the  experiment  coexists  with, 
and  is  functionally  dependent  upon,  the  boundary  layer  build-up  during  the  start-up 
period  of  the  nozzles.  The  boundary  layer  and  the  recompression  shock  system 


- 19  - 


cannot  be  modeled  by  the  inviscid  computational  scheme.  With  the  exception  of  the 
6°  half-nozzle,  the  computations  predict  only  one  recompression  shock  which  gen¬ 
erally  does  not  coincide  with  the  experimental  system,  and  moves  much  slower  than 
the  experimental  shocks  at  low  driver  pressures,  and  much  faster  at  high  driver 
pressures.  In  between,  there  exists  one  pressure  level  at  which  the  computational 
recompression  shock  happens  to  coincide  with  one  of  the  experimental  shocks.  This 
pressure  level  is  different  for  each  divergent  nozzle  angle. 

The  experimental  evidence  suggests  that  the  system  of  multiple  X-type 
recompression  shocks  does  not  move  nearly  as  far  downstream  into  the  expansion 
tunnel  as  does  the  computational  single  recompression  shock  and  decays  faster  than 
the  computation  predicts.  The  computation  also  cannot  predict  the  standing 
recompression  shocks  in  the  throat  channel.  It  predicts  a  moving  shock  instead 
which  weakens  and  dissipates  as  it  approaches  the  exit  of  the  throat  channel.  For 
the  90°  half-nozzle,  the  computation  produces  a  standing  recompression  shock  at  the 
exit  of  the  throat  channel  at  all  pressure  levels  in  addition  to  the  aforementioned 
shock. 


6.  CONCLUSIONS 

Four  convergent-divergent  nozzle  configurations  were  successfully  tested  at 
seven  driver  pressure  ratios  to  study  their  flow  start-up  processes.  However,  the 
field  of  view  for  the  6°  divergent  half-nozzle  and  the  90°  half-nozzle  downstream 
from  the  nozzle  exit  was  too  limited,  and  additional  optical  test  records  are  desirable 
for  a  comparative  study  of  the  flow  start-up  processes  downstream  from  the  exits  of 
the  90°  and  6°  half-nozzles.  The  results  for  the  6°  divergent  half-nozzle  were  not 
discussed  here  because  they  were  too  incomplete  to  contribute  to  the  discussion. 

The  pressure  histories  and  the  optical  flow  records  reveal  that  the  flow  is  very 
sensitive  to  local  disturbances,  e.g.,  flow  asymmetries  caused  by  the  breaking  of  the 
diaphragm,  the  passage  of  the  diaphragm  fragments  through  the  flow  field,  and  the 
boundary-layer  build-up  during  the  start-up  period  of  the  nozzles.  Strong  reflec¬ 
tions  of  the  incident  shock  at  the  exit  corner  of  the  45°  and  90°  half-nozzles  interfere 
with  the  flow  during  the  start-up  process.  A  system  of  multiple  recompression 
shocks  is  generated  by  the  flow  in  the  throat  channel  and  in  the  divergent  section  of 
the  half-nozzle  in  response  to  the  boudary  layer  build-up  along  the  walls  of  the  flow 
channel. 

From  the  observations  and  the  analysis  of  the  available  data,  the  following  con¬ 
clusions  could  be  drawn. 

a)  The  90°  nozzle  effectively  forms  a  divergent  section  by  establishing  its  particular 
flow  expansion  angle  in  response  to  the  particular  driver  pressure  ratio,  but  the 
penalty  to  be  paid  is  a  loss  in  shock  strength  of  about  10  percent  compared  to  the 
shock  recovery  in  the  16°  divergent  half-nozzle  at  the  highest  driver  pressure  ratio. 
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b)  The  shock  formation  period  ranges  from  440±  50p.s  at  the  low  end  of  the  pres¬ 
sure  scale  to  250±  25p.s  at  full  scale  for  all  nozzles.  The  shock  formation  distance  is 
practically  constant  over  the  entire  pressure  range  with  a  mean  value  of  157  mm  and 
a  standard  deviation  of  the  recorded  values  of  o„_  i  =  16  mm. 

c)  The  start-up  period  of  the  45°  nozzle  ranges  from  2.55  ms  at  the  lowest  pressure 
level  (P41  =  4)  to  1.12ms  at  the  highest  pressure  level  (P41  =  188).  The  information 
on  the  start-up  period  for  the  other  nozzles  is  incomplete. 

d)  No  pressure  spikes  at  the  incident  shock  front  of  the  pressure  histories  were 
observed  in  these  experiments  such  as  were  found  in  the  experimental  CEG-LBS 
pressure  histories  or  in  the  BRL-Q1D  computations.  This  is  thought  to  be  due  to 
the  low  area  expansion  ratio  (1 : 7.5)  between  the  throat  and  the  expansion  tunnel. 

e)  The  BRL-Q1D  code  proved  to  be  an  inappropriate  tool  for  studying  the  detailed 
flow  phenomena  of  the  start-up  process  in  nozzles.  It  remains  a  useful  engineering 
tool  for  generating  Oth-order  estimates  of  such  flow  phenomena  and  for  parametric 
design  studies  in  shock  tubes  and  tunnels. 

Based  on  the  present  findings,  it  is  recommended  that  the  experiment  be 
repeated  for  higher  throat/test-section  area  expansion  ratios  of  1:15  and  1:30  with 
redesigned  divergent  nozzle  models  of  16°,  30°,  45°,  and  90°.  The  influence  of  a 
partial  area  discontinuity  at  the  nozzle  exit  should  be  studied  in  this  future  phase  of 
the  experiment.  Since  the  half-nozzle  configuration  did  not  yield  the  desired  one- 
half  symmetrical  flow  pattern  above  the  center  axis  (p.  12),  future  tests  should 
include  at  least  one  full-nozzle  configuration. 

The  90°  half-nozzle  is  of  special  interest  to  the  project  engineer  because  the 
elimination  of  the  divergent  nozzles  from  the  LB/TS  design  saves  a  great  deal  of 
material  cost  in  nozzle  construction.  It  was  found  that  the  flow  through  the  90° 
half-nozzle  with  its  area  discontinuity  at  the  throat  channel  exit  behaves  like  the  flow 
through  the  45°  half-nozzle  but  yields  a  lower  shock  pressure  than  the  latter.  The 
money  saved  by  not  having  an  efficient  nozzle  design  may  be  lost  several  times  over 
by  requiring  a  higher  design  pressure  for  the  drivers  and,  thereby,  thicker  walls  and 
a  stronger  pressure  delivery  system.  Therefore,  45°  divergent  nozzles  are  recom¬ 
mended  for  the  full-scale  LB/TS  to  avoid  additional  losses  in  shock  pressure. 


REFERENCES 

Amann,  H.-O.  "Vorgaenge  beim  Start  einer  ebenen  Reflexionsduese"  (Processes 
during  the  Start-up  of  a  Planar  Reflection  Nozzle).  Report  Nr.  9/68,  Fraunhofer 
Gesellschaft,  Ernst-Mach-Institut,  Freiburg  im  Breisgau,  (West)  Germany,  1968. 

Amann,  H.  O.  "Das  Spike  Problem  bei  der  Simulation  von  Blastwellen"  (The  Spike 
Problem  in  Simulated  Blast  Waves).  EMI-Report  E 13/85,  Fraunhofer  Gesellschaft, 
Ernst-Mach-Institut,  Freiburg  im  Breisgau,  (West)  Germany,  May  1985. 


-  21  - 


Beam,  R.  M.,  and  R.  F.  Warming.  "An  Implicit  Factored  Scheme  for  the  Compres¬ 
sible  Navier-Stokes  Equations  II:/  The  Numerical  ODE  Connection."  Paper  No. 
79-1446,  AIAA  4th  Computational  Fluid  Dynamics  Conference,  Williamsburgh, 
VA,  23-24  July  1979. 

Cadet,  A.,  and  J.  B.  G.  Monzac.  "The  Large-Scale  Nuclear-Blast  Simulator  of  the 
Gramat  Research  Center:  Description  and  Operational  Utilization."  Proceedings  of 
the  Seventh  International  Symposium  on  the  Military  Application  of  Blast  Simula¬ 
tion  (MABS-7),  vol.  1,  paper  no.  1.2,  Medicine  Hat,  Alberta,  Canada,  13-17  July 
1981. 

Coulter,  G.  A.  "Blast  Parametric  Study  Using  a  1:57  Scale  Single  Driver  Model  of  a 
Large  Blast  Simulator",  BRL-MR-3597,  U.  S.  Army  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  MD,  June  1987. 

Emmons,  H.  W.  (ed.)  "Fundamentals  of  Gas  Dynamics."  High  Speed  Aerodynam¬ 
ics  and  Jet  Propulsion,  vol.  Ill,  New  Jersey:  Princeton  University  Press,  1958. 

Ethridge,  N.  H.,  R.  E.  Lottero,  J.  D.  Wortman,  and  B.  P.  Bertrand.  "Computa¬ 
tional  and  Experimental  Studies  of  Blockage  Effects  in  a  Blast  Simulator." 
ARBRL-TR-02564,  U.S.  Army  Armament  Research  and  Development  Center, 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  June  1984. 

Gion,  E,  J.  "A  Multiple  Driver  Shock  Tube  Model  of  a  Large  Blast  Simulator." 
BRL-MR-3757,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving 
Ground,  MD,  May  1989. 

Gion,  E.  J.  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
MD,  Private  Communications,  1989. 

Gratias,  S.,  and  J.  B.  G.  Monzac.  "The  Large-Scale  Nuclear-Blast  Simulator  of  the 
Gramat  Research  Center:  Concept,  Research,  Performance."  Proceedings  of  the 
Seventh  International  Symposium  on  the  Military  Application  of  Blast  Simulation 
(MABS-7),  vol.  1,  paper  no.  1.1,  Medicine  Hat,  Alberta,  Canada,  13-17  July  1981. 

Hisley,  D.  M.,  and  G.  A.  Molvik.  "Axisymmetric  Calculations  for  the  Large 
Blast/rhermal  Simulator  (LEA’S)  Shock  Tube  Configuration."  BRL-TR-2935,  U.S. 
Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  September 
1988. 

Mark,  A.,  K.  O.  Opalka,  C.  W.  Kitchens,  G.  A.  Coulter,  G.  Bulmash  and  C.  N. 
Kingery.  "Simulation  of  Nuclear  Blasts  with  Large-Scale  Shock  Tubes."  Proceedings 
of  the  Eighth  International  Symposium  on  Military  Applications  of  Blast  Simulation 
(MABS-8),  Spiez,  Switzerland,  20-24  June  1983. 

Opalka,  K.  O.  "Large  Blast  and  Thermal  Simulator  Advanced  Concept  Driver 
Design  by  Computational  Fluid  Dynamics",  BRL-TR-3026,  U.S.  Army  Ballistic 


-  22  - 


Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  August  1989. 

Opalka,  K.  O.,  and  R.  J.  Pearson.  "Advanced  Concept  Large  Blast/Thermal  Simu¬ 
lator  CFD  Design  Studies."  Computers  in  Engineering  1989  Conference  Proceed¬ 
ings,  pp.  189-195,  ASME,  New  York,  NY,  August  1989. 

Opalka,  K.  0.,  and  R.  J.  Pearson.  "Real  Time  Flow  Control  in  Large 
Blast/Thermal  Simulators."  Proceedings  of  the  1988  Army  Science  Conference,  Ft. 
Monroe,  Hampton,  VA,  25-28  October  1988. 

Opalka,  K.  O.  "Large  Blast-Wave  Simulators  (LBS)  With  Cold-Gas  Drivers:  Com¬ 
putational  Design  Studies",  BRL-TR-2786,  U.S.  Army  Ballistic  Research  Labora¬ 
tory,  Aberdeen  Proving  Ground,  MD,  March  1987. 

Opalka,  K.  O.,  and  A.  Mark.  "The  BRL-Q1D  Code:  A  Tool  for  the  Numerical 
Simulation  of  Flows  in  Shock  Tubes  with  Variable  Cross-Sectional  Areas,"  BRL- 
TR-2763,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
MD,  October  1986. 

Pearson,  R.  J.  "Large  Blast/Thermal  Simulation."  SAE  Technical  Paper  #  871746, 
Society  of  Automotive  Engineers,  400  Commonwealth  Drive,  Warrendale,  Pennsyl¬ 
vania,  October  1987. 

Pearson,  R.  J.,  K.  O.  Opalka,  and  D.  M.  Hisley.  "Design  Studies  of  Drivers  for  the 
US  Large  Blast/Thermal  Simulator,"  Proceedings  of  the  Ninth  International  Sympo¬ 
sium  on  Military  Applications  of  Blast  Simulation  (MABS-9),  vol.  1,  paper  no.  1.6, 
Southend-on-Sea,  Essex,  England  SS3  9XE,  September  1985. 

Precision  Visuals,  Inc.  "Picsure  User’s  Guide."  6260  Lookout  Road,  Boulder,  CO, 
December  1986. 

Reichenbach,  H.,  and  K.  O.  Opalka.  "An  Optical  Study  of  the  Flow  Start-up  Pro¬ 
cess  in  Four  Convergent-Divergent  Nozzles."  EMI-Report  E  3/90,  Fraunhofer 
Gesellschaft,  Ernst-Mach-Institut,  Freiburg  im  Breisgau,  (West)  Germany,  March 
1990. 

Schardin,  H.  "Das  Toeplersche  Schlieren  Verfahren"  (The  Toepler  Schlieren 
Method).  V.D.I.  Forschungsheft  366,  1934. 

Schraml,  S.  J.,  and  R.  J.  Pearson.  "Small  Scale  Shock  Tube  Experiments  Using  a 
Computer  Controlled  Active  Rarefaction  Wave  Eliminator."  BRL-TR-3149,  U.S. 
Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  September 
1990. 


-23  - 


(intentionally  left  blank) 


-24- 


-  25  - 


LEGEND: 


-Active  RWE  Boundary 
■Driver  Pressure  Boundary 


.600  kT 


Positive-Phase  Duration,  s 


(a)  Based  on  Positive  Phase  Duration 


h  fcT/V  / 

/  /  /■ 

/  /10  kT 


/  / 
/  / 


'100  kT 


600  kT 


LEGEND: 

- Active  RWE  Boundary 


Driver  Pressure  Boundary 


Dynamic  Pressure  Impulse,  kPa-s 


(b)  Based  on  Dynamic  Pressure  Impulse 


Figure  2: 


sign  hnve 


12  9 


LEGEND 
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Figure  3:  Large  Blast/Thermal  Simulator  Design  Concept. 
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Figure  6:  EMI  Planar  Driver  and  Nozzle  Design. 


(b)  Diaphragm  Section 
Figure  6:  (concluded) 
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(d)  6°  Half-Nozzle 
Figure  7:  (concluded) 
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Figure  9:  Triggering  Circuit  (EMI 
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a)  Shadowgraph  Sequence  in  Early  Phase  of  Flow  Start-up 


b)  Illustration  of  Major  Flow  Phenomena  in  Early  Phase 
Figure  10:  (continued) 
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c)  Shadowgraph  Sequence  in  L 
Figure  10:  (< 
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Phase  of  Flow  Start-up 


d)  Illustration  of  Major  Flow  Phenomena  in  Late  Phase 
Figure  10:  (concluded) 
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a)  Driver  Pressure  Ratio  P41  =  4 

Figure  11:  Two  Sequences  of  Selected  Shadowgraphs  for  the  90°  Half-Nozzle. 
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b)  Driver  Pressure  Ratio  P41  =  188 
Figure  11:  (concluded) 
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b)  Short  Time  Seale 

Figure  15:  Static  Overpressure  History  for  the  16°  Half-Nozzle  at  Driver  Pressure 
Ratio  =  14.  860mm  Downstream  from  the  Diaphragm. 
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a)  Shock  Strength  vs.  Driver  Pressure  Ratio 


b)  Shock  Strength  vs.  Shock  Mach  Number 
Figure  19:  Shock  Strength  at  Gage  #  37  (^Omm  from  the  Diaphragm). 
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Shock  Strength,  PI(  Shock  Strength,  P; 


a)  Shock  Strength  vs.  Driver  Pressure  Ratio 


b)  Shock  Strength  vs.  Shock  Mach  Number 
Figure  20:  Shock  Strength  at  Gage  #  86  (860mm  from  the  Diaphragm"). 
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Shock  Strength,  PI(  Shock  Strength, 


a. 


a)  Shock  Strength  vs.  Driver  Pressure  Ratio 


b)  Shock  Strength  vs.  Shock  Mach  Number 
Figure  21:  Shock  Strength  at  Gage  #  99  (989mm  from  the  Diaphragm). 
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Figure  A-2.  6°  Nozzle  at  P4i=  14 
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Figure  A-3:  6°  Nozzle  at  P4i=38. 
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Figure  A-5.  6°  Nozzle  at  P 
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Figure  A-6.  6°  Nozzle  at  P< 
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Figure  A-7.  6°  Nozzle  at  P4  L88 
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Figure  A-8.  16°  Nozzle  at  P4i=  4 
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Figure  A-9.  16°  Nozzle  at  P«=  14 
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Figure  A- 10.  16°  Nozzle  at  P. 
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Figure  A— 11.  16®  Nozzle  at  P. 
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Figure  A- 12.  16°  Nozzle  at  P 
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Figure  A— 13.  16°  Nozzle  at  P 
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Figure  A- 14.  16°  Nozzle  at  P4i=  188 
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Figure  A- 15:  45°  Nozzle  at  P«=  4. 
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Figure  A- 16.  45°  Nozzle  at  P4i=  14 
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Figure  A- 17:  45°  Nozzle  at  P4i=38. 
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Figure  A- 18.  45°  Nozzle  at  P„i=  55 


-72- 


Shot  No. 1 3308  (0-200mm) 


Shot  No. 13319  (100-300mm) 


— Incident  Shock 
+  Contact  Surface 

★  1st  Recompression  Shock 
■  2nd  Recompression  Shock 

-  o-  First  Diaphragm  Fragment 

-  o-  Last  Diaphragm  Fragment 

a  2nd  Shock  in  Throat  Ch’l  o^ 

♦  X-Type  Throat  Shock  . 


/  © 
o'  n 


Distance  (mm) 


Figure  A- 19.  45°  Nozzle  at  P 
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Figure  A-20.  45°  Nozzle  at  P. 
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Figure  A-21.  45°  Nozzle  at  P4i=  188 
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Figure  A-22.  90°  Nozzle  at  P4i~  4 
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Figure  A-23.  90*  Nozzle  at  P4i=  14 
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Figure  A— 26:  90°  Nozzle  at  P.ii=80. 
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Figure  A-27:  90°  Nozzle  at  P 
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Figure  A— 28:  90°  Nozzle  at  P, 


APPENDIX  B:  X-t  DIAGRAMS  WITH  Q1D  COMPARISONS 


-  83  - 


(intentionally  left  blank) 


-  84  - 
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Figure  B-2:  Case  N06/014 
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Figure  B-3:  Case  N06/038 
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Figure  B-4:  Case  N06/080 
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Figure  B-6:  Case  N16/GJ.4. 
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Figure  B-7:  Case  N16/038 
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Figure  B-8:  Case  N 16/080. 
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Figure  B-9:  Case  N45/004. 
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Figure  B-10:  Case  N45/014 
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Figure  B-ll:  Case  N45/038 


Time 


-  96  - 


Time  (ms) 


Shots  No.  13249  &  13254  (0-200mm) 


•  •  Incident  Shock 

•  -  Q1D  Incident  Shock 

0  +  Contact  Surface 

- Q1D  Contact  Surface 

*  1st  Recompression  Shock 

*  - Q1D  Recompression  Shock 

•  •  2nd  Recompression  Shock 

■  3rd  Recompression  Shock 

*  x  3rd  Shock  in  Throat  Ch’l 

'  a  2nd  Shock  in  Throat  Ch’^ ' 
- Q1D  Shock  in  Throat  CJ?'I 

♦  Oblique  Throat  Shock  f 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


Distance  (mm) 


Figure  B-13:  Case  N90/004 
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Figure  B— 14:  Case  N90/014 
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Figure  B— 15:  Case  N90/038. 
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Figure  B-16:  Case  N90/080. 
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SMCCR-MSI 

Dir,  VLAMO 
ATTN:  AMSLC-VL-D 


1  Director  10  Dir,  BRL 

U.S.  Army  Aviation  Research  ATTN:  SLCBR-DD-T 

and  Technology  Activity 
ATTN:  SAVRT-R  (Library) 

M/S  219-3 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 
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1  Director  of  Defense  Research  & 
Engineering 
ATTN:  DD/TWP 
Washington,  DC  20301 

1  Assistant  Secretary  of  Defense 
(Atomic  Energy) 

ATTN:  Document  Control 
Washington,  DC  20301 

1  Chairman 

Joint  Chiefs  of  Staff 
ATTN:  J-5,  R&D  Division 
Washington,  DC  2C301 

2  Deputy  Chief  of  Staff  for  Operations 

and  Plans 

ATTN:  Technical  Library 

Director  of  Chemical 
and  Nuclear  Operations 
Department  of  the  Army 
Washington,  DC  20310 

1  European  Research  Office 

USARDSG  (UK) 

ATTN:  Dr.  R.  Reichenbach 
Box  65 

FPO  New  York  09510-1500 

1  Director 

Defense  Advanced  Research  Projects 
Agency 

ATTN:  Tech  Lib 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

2  Director 

Federal  Emergency  Management 
Agency 

ATTN:  Public  Relations  Office 
■  Technical  Library 
Washington,  DC  20472 

1  Chairman 

DOD  Explosives  Safety  Board 
Room  856-C 
Hoffman  Bldg.  1 
2461  Eisenhower  Avenue 
Alexandria,  VA  22331-0600 


No.  of 

Copies  Organization 
1  Director 

Defense  Intelligence  Agency 
ATTN:  DT-2/Wpns  &  Sys  Div 
Washington,  DC  20301 

1  Director 

National  Security  Agency 
ATTN:  R15,  E.  F.  Butala 
Ft.  George  G.  Meade,  MD  20755 

7  Director 

Defense  Nuclear  Agency 
ATTN:  CSTI,  Tech  Lib 
DDIR 

DFSP,  Ullrich 

NANS 

OPNA 

SPSD 

SPTD 

Washington,  DC  20305 

3  Commander 

Field  Command,  DNA 
ATTN:  FCPR 

FCTMOF 

NMHE,  CDR  Lund 
Kirtland  AFB,  NM  87115 

10  Centra!  Intelligence  Agency 
DIR/DB/Standard 
ATTN:  GE-47HQ 
Washington,  DC  20505 

5  Director 

US  Army  Harry  Diamond  Laboratories 
ATTN:  SLCHD-NW-RA, 

L.  Belliveau 
SLCHD-NW-P, 

Patnaik 

Abe 

Corrigan 

SLCHD-TA-L,  Tech  Lib 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 

2  Commander,  USACECOM 
ATTN:  AMSEL-RD 

AMSEL-RO-TPPO-P 
Fort  Monmouth,  NJ  07703-5301 
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1  Commander,  USACECOM 
R&D  Technical  Library 
ATTN:  ASQNC-ELC-I-T,  Myer  Center 
Fort  Monmouth,  NJ  07703-5301 

1  Director 

US  Army  Missile  and  Space 
Intelligence  Center 
ATTN:  AIAMS-YDL 
Redstone  Arsenal,  AL  35898-5500 

1  Commander 

US  Army  Foreign  Science  and 
Technology  Center 
ATTN:  Research  &  Data  Branch 
220  7th  Street,  NE 
Charlottesville,  VA  22901 

1  Director 

US  Army  TRAC  -  Ft.  Lee 
ATTN:  ATRC-L,  R.  Cameron 
Fort  Lee,  VA  23801-6140 

1  Commander 

US  Army  Materials  Technology 
Laboratory 

ATTN:  AMXMR-ATL 
Watertown,  MA  02172-0001 

2  Commander 

US  Army  Strategic  Defense 
Command 

ATTN:  CSSD-H-MPL,  Tech  Lib 
CSSD-H-XM,  Dr.  Davies 
PO  Box  1500 
Huntsville,  AL  35807 

2  Commander 

US  Army  Natick  Research  and 
Development  Center 
ATTN:  AMDNA-D,  Dr.  D.  Sieling 
STRNC-UE,  J.  Calligeros 
Natick,  MA  01762 

1  Commander 

US  Army  Engineer  Division 
ATTN:  HNDED-FD 
PO  Box  1500 
Huntsville,  AL  35807 


No.  of 

Copies  Organization 

3  Commander 

US  Army  Corps  of  Engineers 
Waterways  Experiment  Station 
ATTN:  CAWES-SS-R,  J.  Watt 
CAWES-SE-R,  J.  Ingram 
CAWES-TL,  Tech  Lib 
PO  Box  631 

Vicksburg,  MS  39180-0631 

1  Commander 

US  Army  Corps  of  Engineers 
Fort  Worth  District 
ATTN:  CESWF-PM-J, 

R.  Timmermins 
P.O.  Box  17300 
Forth  Worth,  TX  76102-0300 

1  Commander 

US  Army  Research  Office 
ATTN:  SLCRO-D 
POBox  12211 

Research  Triangle  Park,  NC  27709-2211 
3  Commander 

US  Army  Nuclear  &  Chemical  Agency 
ATTN:  ACTA-NAW 
MONA-WE 
Tech  Lib 

7500  Backlick  Rd.,  Bldg.  2073 
Springfield,  VA  22150 

1  Director 

HQ,  TRAC  RPD 

ATTN:  ATRC-RPR,  Mr.  Radda 

Fort  Monroe,  VA  23651-5143 

1  Director 

TRAC-WSMR 

ATTN:  ATRC-WC,  Mr.  Kirby 
White  Sands  Missile  Range,  NM 
88002-5502 

1  Director 
TRAC-FLVN 
ATTN:  ATRC 

Fort  Leavenworth,  KS  66027-5200 
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Copies  Organization 


1  Commander 

US  Army  Test  &  Evaluation 
Command 

Nuclear  Effects  Laboratory 
ATTN:  STEWS-TE-NO, 

Dr.  J.  L.  Meason 
PO  Box  477 

White  Sands  Missile  Range,  NM 
88002 

2  Chief  of  Naval  Operations 
ATTN:  OP-03EG 

OP-985F 

Department  of  the  Navy 
Washington,  DC  20350 

1  Director 

Strategic  Systems  Projects  Office 
ATTN:  NSP-43,  Tech  Library 
Department  of  the  Navy 
Washington,  DC  20360 

1  Commander 

Naval  Electronic  Systems  Command 
ATTN:  PME  117-21 A 
Washington,  DC  20360 

1  Commander 

Naval  Facilities  Engineering  Command 
ATTN:  Technical  Library 
Washington,  DC  20360 

1  Commander 

Naval  Sea  Systems  Command 
ATTN:  Code  SEA-62R 
Department  of  the  Navy 
Washington,  DC  20362-5101 

2  Office  of  Naval  Research 
ATTN:  Dr.  A.  Faulstick,  Code  23 
800  N.  Quincy  Street 
Arlington,  VA  22217 

1  Officer-in-Charge 

Naval  Construction  Battalion  Center 
Civil  Engineering  Laboratory 
ATTN:  Tech  Lib,  Code  L06C/L08A 
Port  Hueneme,  CA  93041 


1  Commanding  Officer 

Naval  Civil  Engineering  Laboratory 
ATTN:  Code  L51,  J.  Tancreto 
Port  Hueneme,  CA  93043-5003 

1  Commander 

Naval  Surface  Warfare  Center 
ATTN:  Code  DX-21 , '  ‘brary 
Dahlgren,  VA  22448-5000 

1  Commander  (Code  533) 

Naval  Weapons  Center 
Tech  Library 

China  Lake,  CA  93555-6001 

1  Commander 

Naval  Weapons  Evaluation  Fac 
ATTN:  Document  Control 
Kirtland  AFB,  NM  87117 

1  Commander 

Naval  Research  Laboratory 
ATTN:  Code  2027,  Tech  Library 
Washington,  DC  20375 

2  Air  Force  Armament  Laboratory 
ATTN:  AFATL/DOIL 

AFATL/DLYV 

Eglin  AFB,  FL  32542-5000 

1  AFESC/RDCS 

ATTN:  Paul  Rosengren 
Tyndall  AFB,  FL  32403 

1  RADC  (EMTLD/Docu  Library) 

Griff iss  AFB,  NY  13441 

3  Air  Force  Weapons  Laboratory 
ATTN:  NTE 

NTED 

NTES 

Kirtland  AFB,  NM  87117-6008 
1  AFIT 

ATTN:  Tech  Lib,  Bldg.  640/B 
Wright-Patterson  AFB,  OH  45433 

1  ALASCF 

ATTN:  J.  Levine 

Edwards  AFB,  CA  93523-5000 
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1  AL/TSTL  (Tech.  Lib) 

ATTN:  J.  Lamb 

Edwards  AFB,  CA  93523-5000 

1  FTD/NIIS 

Wright-Patterson  AFB,  OH  45433 


Director 

Los  Alamos  National  Laboratory 
ATTN:  Mr.  Th.  Dowler,  MS-F602 

Doc  Control  for  Reports  Library 
PO  Box  1663 
Los  Alamos,  NM  87545 


1  US  Department  of  Energy 
Idaho  Operations  Office 

ATTN:  Spec  Programs,  J.  Patton 
785  DOE  Place 
Idaho  Falls,  ID  83402 

2  Director 

Idaho  National  Engineering  Laboratory 
EG&G  Idaho  Inc. 

ATTN:  Mr.  R.  Guenzler,  MS-3505 
Mr.  R.  Holman,  MS-3510 
PO  Box  1625 
Idaho  Falls,  ID  83415 

1  Commander 

David  Taylor  Research  Center 
ATTN:  Code  522,  Tech  Info  Ctr 
Bethesda,  MD  20084-5000 

1  Officer  in  Charge 

White  Oak  Warfare  Center 
Detachment 

ATTN:  Code  E232,  Tech  Library 
10901  New  Hampshire  Avenue 
Silver  Spring,  MD  20903-5000 

1  Commanding  Officer 

White  Oak  Warfare  Center 
ATTN:  Code  WA501,  NNPO 
Silver  Spring,  MD  20902-5000 

1  Director 

Lawrence  Livermore  National  Laboratory 
ATTN:  Tech  Info  Dept  L-3 
PO  Box  808 
Livermore,  CA  9455Q 

1  Director 

NASA-Langiey  Research  Center 
ATTN:  Tech  Lib 
Hampton,  VA  23665 


3  Director 

Sandia  National  Laboratories 
ATTN:  Doc  Control  3141 

Mr.  C.  Cameron,  Div  6215 
Mr.  A.  Chabai,  Div  7112 
PO  Box  5800 

Albuquerque,  NM  87185-5800 
1  Director 

Sandia  National  Laboratories 

Livermore  Laboratory 

ATTN:  Doc  Control  for  Tech  Library 

PO  Box  969 

Livermore,  CA  94550 

1  Director 

National  Aeronautics  and  Space 
Administration 

ATTN:  Scientific  &  Tech  Info  Fac 
PO  Box  8757,  BWI  Airport 
Baltimore,  MD  21240 

1  Director 

NAS  A- Ames  Research  Center 
Applied  Computational 
Aerodynamics  Branch 
ATTN:  Dr.  T.  Holtz,  MS  202-14 
Moffett  Field,  CA  94035 

2  Applied  Research  Associates,  Inc. 
ATTN:  N.  H.  Ethridge 

J.  Keefer 
PO  Box  548 
Aberdeen,  MD  21001 

1  Aerospace  Corporation 
ATTN:  Tech  Info  Services 
PO  Box  92957 
Los  Angeles,  CA  90009 
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1  Agbabian  Associates 
ATTN:  M.  Agbabian 
250  North  Nash  Street 
El  Segundo,  CA  90245 

1  Applied  Research  Associates,  Inc. 

ATTN:  R.  L.  Guice 
71 14  West  Jefferson  Ave.,  Suite  305 
Lakewood,  CO  80235 

1  Black  &  Veatch,  Engineers-Architects 

ATTN:  H.  D.  Laverentz 
1500  Meadow  Lake  Parkway 
Kansas  City,  MO  64114 

1  The  Boeing  Company 
ATTN:  Aerospace  Library 
PO  Box  3707 
Seattle,  WA  98124 

1  California  Research  &  Technology,  Inc. 

ATTN:  M.  Rosenblatt 
20943  Devonshire  Street 
Chatsworth,  CA  91311 

1  Carpenter  Research  Corporation 

ATTN:  H.  Jerry  Carpenter 
27520  Hawthorne  Blvd.,  Suite  263 
PO  Box  2490 

Rolling  Hills  Estates,  CA  90274 

1  Dynamics  Technology,  Inc. 

ATTN:  D.  T.  Hove 

21311  Hawthorne  Blvd.,  Suite  300 

Torrance,  CA  90503 

1  EATON  Corporation 

Defense  Valve  &  Actuator  Div. 

ATTN:  Mr.  J.  Wada 
2338  Alaska  Ave. 

El  Segundo,  CA  90245-4896 

1  Goodyear  Aerospace  Corporation 

ATTN:  R.  M.  Brown,  Bldg.  1 
Shelter  Engineering 
Litchfield  Park,  A Z  85340 


No.  of 

Copies  Organization 

1  Altiant  Techsystems,  Inc. 

ATTN:  Roger  A.  Rausch  (MN48-3700) 
7225  Northland  Drive 
Brooklyn  Park,  MN  55428 

4  Kaman  AviDyne 

ATTN:  Dr.  R.  Ruetenik  (2  cys) 

Mr.  S.  Criscione 
Mr.  R.  Milligan 
83  Second  Avenue 
Northwest  Industrial  Park 
Burlington,  MA  01803 

3  Kaman  Sciences  Corporation 
ATTN:  Library 

P.  A.  Ellis 
F.  H.  Shelton 
PO  Box  7463 

Colorado  Springs,  CO  80933-7463 

1  Kaman  Sciences  Corporation 
ATTN:  Mr.  F.  W.  Balicki 

6400  Uptown  Boulevard  N.E.,  Suite  300 
Albuquerque,  NM  871 1 0 

2  Kaman-TEMPO 
ATTN:  DASIAC 

Don  Sachs 
PO  Drawer  1479 
816  State  Street 
Santa  Barbara,  CA  93102-1479 

1  Lockheed  Missiles  &  Space  Co. 

ATTN:  J.  J.  Murphy 

Dept.  81-11,  Bldg.  154 
PO  Box  504 
Sunnyvale,  CA  94086 

2  McDonnell  Douglas  Astronautics 

Corporation 

ATTN:  Robert  W.  Halprin 
K.A.  Heinly 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 ' 

2  The  Ralph  M.  Parsons  Company 

ATTN:  Tyler  M.  Jackson, 

LB/TS  Project  Manager 
100  West  Walnut  Street 
Pasadena,  CA  91124 


No.  of 
Copies 


1 


3 


1 


1 


1 


1 


1 


Organization 


No.  of 

Copies  Organization 


R&D  Associates 
ATTN:  Technical  Library 
Dr.  Allan  Kuhl 
PO  Box  9695 

Marina  Del  Rey,  CA  90291 

R&D  Associates 
ATTN:  G.  P.Ganong 
PO  Box  9377 
Albuquerque,  NM  87119 

Science  Applications,  Inc. 

ATTN:  Division  164,  M  ST-302 
W.  Layson 
John  Cockayne 
PO  Box  1303 
1710  Goodridge  Drive 
McLean,  VA  22102 

Science  Applications  International 
Corp. 

ATTN:  Mr.  J.  Guest 
2301  Yale  Blvd.  SE,  Suite  E 
Albuquerque,  NM  87106 

Science  Applications  International 
Corporation 
ATTN:  Mr.  N.  Sinha 
501  Office  Center  Drive,  Apt.  420 
Ft.  Washington,  PA  19034-3211 

Sparta,  Inc. 

Los  Angeles  Operations 
ATTN:  I.  B.  Osofsky 
3440  Carson  Street 
Torrance,  CA  90503 

Sverdrup  Technology,  Inc. 

ATTN:  R.  F.  Starr 
PO  Box  884 
Tullahoma,  TN  37388 


2  S-CUBED 

A  Division  of  Maxwell 
Laboratories,  Inc. 

ATTN:  C.  E.  Needham 
Dr.  Lynn  Kennedy 
2501  Yale  Blvd.,  SE 
Albuquerque,  NM  87106 

3  S-CUBED 

A  Division  of  Maxwell 
Laboratories,  Inc. 

ATTN:  Technical  Library 
R.  Duff 
K.  Pyatt 
PO  Box  1620 
La  Jolla,  CA  92037-1620 

1  Texas  Engineering  Experiment 
Station 

ATTN:  Dr.  D.  Anderson 

301  Engineering  Research  Center 

College  Station,  TX  77843 

1  Thermal  Science,  Inc. 

ATTN:  R.  Feldman 
2200  Cassens  Dr. 

St.  Louis,  MO  63026 

1  Battelle 

ATTN:  TACTEC  Library,  J.  N.  Higgins 
505  King  Avenue 
Columbus,  OH  43201-2693 

1  California  Institute  of  Technology 

ATTN:  T.  J.  Ahrens 
1201  E.  California  Blvd. 

Pasadena,  Ca  91109 

1  Massachusetts  Institute  of  Technology 

ATTN:  Technical  Library 
Cambridge,  MA  02139 


SRI  International 
ATTN:  Dr.  G.  R.  Abrahamson 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 


Massachusetts  Institute  of  Technology 
Aeroelastic  and  Structures  Research 
Laboratory 

ATTN:  Dr.  E.  A.  Witmer 
Cambridge,  MA  02139 
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2  New  Mexico  Engineering  Research 

Institute  (CERF) 

University  of  New  Mexico 
ATTN:  Dr.  J.  Leigh 
Dr.  R.  Newell 
PO  Box  25 

Albuquerque,  NM  87131 

1  Northrop  University 

ATTN:  Dr.  F.  B.  Safford 
5800  W.  Arbor  Vitae  St. 

Los  Angeles,  CA  90045 

3  Southwest  Research  Institute 
ATTN:  C.  Anderson 

S.  Mullin 
A.  B.  Wenzel 
PO  Drawer  28255 
San  Antonio,  TX  78228-0255 

1  State  University  of  New  York 
Mechanical  &  Aerospace  Engineering 
ATTN:  Dr.  Peyman  Givi 

Buffalo,  NY  14260 

2  Denver  Research  Institute 
ATTN:  Mr.  J.  Wisotski 

Technical  Library 
P.O.  Box  10758 
Denver,  CO  80210 

1  Stanford  University 

ATTN:  Dr.  D.  Bershader 
Durand  Laboratory 
Stanford,  CA  94305 

1  AEDC 

ATTN:  R.  McAmis,  Mail  Stop  980 
Arnold  AFB.TN  37389 


Aberdeen  Proving  Ground 

1  Cdr,  USATECOM 

ATTN:  AMSTE-TE-F,  L.  Teletski 

1  Cdr,  USATHMA 
ATTN:  AMXTH-TE 

1  Cdr,  USACSTA 
ATTN:  STECS-LI 
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3  Director 

Fraunhofer-lnstitut  fOr  Kurzzeitdynamik 
Ernst-Mach-lnstitut 
Eckerstrasse  4 
D-7800  Freiburg  (Br.) 

Germany 


(Intentionally  left  blank) 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number  BRL-TR-3215 _ Date  of  Report  March  1991 _ 

2.  Date  Report  Received .  . . 


3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 
for  which  the  report  will  be  used.) _  _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ _ _ — 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.)  _  _ 


Name 

CURRENT  Organization 

ADDRESS _ 

Address 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 

OLD  Organization 

ADDRESS  _ 

Address 

City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


FOLD  HERE 


Department  of  the  Army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21  Co  -5066 

OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  No  0001.  APG,  MD 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T' 

Aberdeen  Proving  Ground,  MD  21005-^)89 


FOLD  HERE 


